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Executive Summary 

This report describes the research performed in the two year period be- 
tween June 1988 and May 1990, under the NASA Johnson Space Center 
Grant NAG 9-275 to Texas A&M University. The research results and 
findings of the two major tasks of the project are included. Namely, the 
report details the design and evaluation of an Attitude Control/Momentum 
Management System for highly asymmetric spacecraft configurations, and 
the preliminary development and application of a nonlinear control sys- 
tem design methodology for tracking control of uncertain systems, such as 
spacecraft payload pointing systems. Control issues relevant to both linear 
and nonlinear rigid-body spacecraft dynamics are addressed, whereas any 
structural flexibilities are not taken into consideration. 

Results from the first task indicate that certain commonly used sim- 
plifications in the equations of motions result in unstable attitude control 
systems, when used for highly asymmetric spacecraft configurations. An ap- 
proach is suggested circumventing this problem. Additionally, even though 
preliminary results from the second task are encouraging, the proposed 
nonlinear control system design method requires further investigation prior 
to its application and use as an effective payload pointing system design 
technique. 
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Part I 


Project Description 


Chapter 1 

Overview of the Project 


In Part II of this report an attitude control/momentum management (ACMM) 
system is proposed for large angle torque-equilibrium- attitude (TEA) con- 
figurations of the Space Station Freedom (SSF). The fully coupled equa- 
tions of motion for the SSF are linearized around an arbitrary TEA and 
the linear-quadratic-regulator (LQR) method with pole placement is used 
to synthesize the full-state feedback gains. The actuators used are control 
moment gyros (CMGs) and both the gyroscopic and the gravity gradient 
torques are utilized in stabilizing the spacecraft, in the presence of exter- 
nal cyclic and secular disturbances. As it is anticipated during the earlier 
flights of the SSF and during some of the Mobile Remote Manipulator Sys- 
tem (MRMS) operations, the contributions of the cross-products of inertia 
on the overall spacecraft attitude determination will be significant, result- 
ing in large pitch angle TEAs. Therefore, a fully coupled control system 
will be required for these earlier flights and for certain modes of operation. 
The proposed gain-scheduled controller can be utilized for attitude control 
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and momentum management of any flight during the assembly sequence 
once the CMGs are available and/or during the MRMS operations. As 
shown in this paper, for certain SSF configurations characterized by large 
pitch angle TEAs, failure to include the contributions of the cross-products 
of the inertia terms in the linearization of the SSF dynamics results in the 
design of unstable controllers. These instabilities can only be observed on 
the nonlinear simulator, with the linearized closed-loop system appearing 
stable. The proposed gain- scheduled controller has been tested with iner- 
tia matrix element variations from Flight ^3 (MB3) to Assembly Complete 
(AC) configurations, corresponding to up to 1600% changes in some cross- 
inertia terms. Computer simulation results indicate that because of the 
slowly changing dynamics of the SSF, a gain-scheduled controller provides 
satisfactory form of adaptation, as compared to other adaptive schemes, 
requiring more complex control architectures. 

Part III of this report is focused on servo- tracking in the context of the 
ability to follow a desired output, within a prespecified error bound, over 
its entire duration; that is, to obtain a closed-loop system which has the 
ability to track a reference input vector while rejecting a class of external 
disturbances, to the extent required by the error bound. This type of 
tracking is referred to as tracking in the sense of spheres. It has been shown 
in the literature that a solution to this tracking problem may be obtained by 
a quantitative pole placement (QPP) approach. In particular, a nonlinear 
observer-based state feedback controller will guarantee the desired tracking 
performance, provided state feedback gains can be found such that a specific 
induced norm of a linear closed-loop operator can be made sufficiently small. 
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In this report a systematic procedure is given for solving the QPP prob- 
lem when the induced norm is the Loo - norm. In particular, a Lure type 
system is considered with the input and output spheres measured with re- 
spect to the Loo - norm, and an algorithm is developed on the basis of 
minimizing the Hoo - norm of a linear closed-loop operator for solving the 
QPP problem. Relevant design criteria axe derived by first embedding the 
problem in the L2 function space, leading to an if 00 - minimization problem 
which is followed by an exponential weighting technique for estimating the 
tracking error bounds in terms of Loo - norms. The approach to this prob- 
lem is based on fixed-point techniques, the multivariable circle criterion, 
and ffoo-optimal control. The controller gains satisfying the design criteria 
and the QPP are obtained by solving a Riccati equation and a Lyapunov 
equation. 

The controller synthesis methodology for servo- tracking in uncertain, 
nonlinear systems of the Lure type presented in the report is applied to a 
simple single input-single output, nonlinear, uncertain system and a two 
degree of freedom robotic manipulator. The manipulator tip is required 
to carry an unknown payload, lying between a given upper and a lower 
bound, along a semi-circular trajectory to a prespecified accuracy. The 
resulting controller resembles a robust computed torque arrangement with 
the feedback gains chosen specifically to provide the necessary reduction 
in sensitivity due to possible variations in the payload. Simulations results 
confirm the effectiveness of the synthesis technique. 


Part II 


Attitude Control/Momentum 
Management of Advanced 
Space Vehicles 



Chapter 2 
Introduction 


The Space Station Freedom (SSF) will be the first permanently manned 
spacecraft to be assembled and operated by NASA. Currently, the schedule 
calls for the First Element Launch (FEL) in March of 1995 and the Assem- 
bly Complete (AC) configuration towards the last quarter of 1999. Two of 
the planned uses of the SSF are as a transportation node for future space 
missions, such as the Lunar and Mars Initiative, and as an experimental 
facility for zero-gravity materials processing and life sciences research. As a 
result of these program objectives, there is a need for the SSF environment 
to be effectively controlled. That is attitude control during assembly, nor- 
mal flight and other operating conditions, such as operation of the Mobile 
Remote Manipulator System (MRMS). 

Even though during the first few flights of the SSF reaction jets will 
be used as actuators, Control-Moment- Gyros (CMGs) will be utilized as 
the primary actuating devices for most of the assembly sequence and sub- 
sequent man-tended operations. Additionally, the gravity gradient torque 
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will be used to unload the spacecraft momentum from the CMGs and bring 
it to its nominal level. Previous studies have presented the trade-offs in se- 
lecting the various momentum management techniques [1]. In the current 
study, the momentum management approach proposed in reference [1] is 
used, however, only the attitude hold mode of the SSF operation is exam- 
ined. The designed control law uses proportional-derivative (PD) for atti- 
tude control and proportional-integral (PI) for momentum management. 

In addressing the attitude hold of the SSF, no explicit design specifi- 
cations are placed upon the control system, other than the desired closed- 
loop system stability. The ultimate objective of this work is the application 
of direct nonlinear control system design techniques to the SSF Attitude 
Control/Momentum Management (ACMM) system. In order, however, to 
determine the benefits, if any, of using a more complex nonlinear controller, 
a gain-scheduled ACMM system is designed for comparison purposes. The 
gain- scheduled control system, a simplified and open-loop form of adaptive 
control, will allow attitude control and momentum management through- 
out the entire operating regime of the SSF dynamics. In order to design 
a gain-scheduled controller, however, it is necessary to have the linearized 
system dynamics at various operating points. These will then be used along 
with linear control synthesis techniques to arrive at the various gain sets 
which will in turn be fitted to obtain an analytic form for the gain matrix 
elements as a function of a system variable, called the scheduling variable. 

This paper addresses the linearization of the equations of motion around 
an arbitrary operating point, which in this case is the torque-equilibrium 
attitude (TEA). It is shown that for some SSF configurations, resulting 
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in large pitch angle TEAs, it is necessary to consider the full coupling of 
the roll/yaw with the pitch equations, through the contributions of the 
cross-products of inertia terms. Computer simulations show that use of 
linearized models without the roll/yaw and pitch coupling in the design of 
the ACMM system, result in closed-loop systems that exhibit instabilities 
for large pitch angle TEA configurations of the SSF. Such configurations 
will be encountered during the assembly sequence and during some of the 
MRMS operations [2]. 

Part II of this report is organized in the following manner. Chapter 3 
presents the complete governing equations of motion for the SSF under the 
assumption of a single rigid body. Following, in chapter 4 is the derivation 
of the linearized equations of motion for an arbitrary TEA resulting from 
possible asymmetric configurations where the control axes are misaligned 
with the body axes. Chapter 5 presents the relevant issues in the design of 
an ACMM control law, as well as the Linear Quadratic Regulator (LQR) 
used in this study for the various SSF configurations. The section concludes 
with the gain-scheduling approach used in this study. Chapter 6 includes 
some of the simulation results necessary to substantiate the major points 
of this study. Chapter 7 summarizes the findings of this study with some 
concluding remarks. 



Chapter 3 

Nonlinear Dynamics of the 
Space-Station Rigid-Body 
Motion 


The SSF is expected to operate in a circular orbit at an altitude of approx- 
imately 200 nm, with a period of approximately 5700 seconds. In reality 
SSF represents a flexible multi-body system, subject to various external 
disturbances which may not be accurately known. During the assembly 
sequence the spacecraft will undergo increasingly complex modifications as 
a results of the attached solar panels, thermal radiators, payloads, pressur- 
ized modules, etc. As a result, there will be an increased interaction among 
the various structural components and the spacecraft guidance, navigation 
and control system. In order to guarantee safe and reliable operation, the 
ACMM system should be able to keep the attitude excursions within some, 
yet unspecified, bounds while using only the available control torque and 
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without saturating the CMGs. Therefore, there is a need to systematically 
design ACMM systems which will have predictable and reliable operation. 

Assuming single rigid-body motion, the nonlinear equations of motion 
for the gravity gradient stabilized SSF, with body-fixed control axes, can 
be expressed as follows [1,3]: 

I«(f) = r gyT0 {t ) + r gg (t) + (w (t) - u(<)) (3-1) 

where w(t) = [u>i (<), w 2 (t), w 3 (t)] r are the body-axis components of the 
absolute angular velocity of the spacecraft, with the subscripts (1,2,3) rep- 
resenting the (roll, pitch, yaw) axes, respectively. The matrix I contains 
the principal moments of inertia, /,,, i = 1,2,3, as well as the products of 
inertia i = 1,2,3 and j = 1,2,3 with i / j. It is also assumed that 
I ti = Iji, i,j = 1,2,3. The first two terms on the right-hand-side of 
equation (3.1) represent the gyroscopic and gravity gradient torques. The 
vectors w(<) = [iOj(t), u> 2 (i), tt> 3 (f)] T and u(<) = [ui(t), u 2 (0> u 3(0] 
represent the body-axis components of the external disturbance torque and 

the CMG control torque, respectively. 

The gyroscopic torque term is expressed in terms of the absolute angular 

velocities as follows: 



o “^3(0 ^2(0 


III I \2 *^13 


wi(t) 

Tgyro(t) — 

U>3 (t) 0 — 


J21 I22 ^23 


CJ2 (t) 




I3I I32 133 


_ ^ 3 (t) 


(3.2) 

The gravity gradient torque term is expressed in terms of the body-axis 
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angles with respect to LVLH as: 



0 “C3(<) C 2 (i) 


ill il2 il3 


1 

r> 

1 

Tgg(t) = 3n 2 

c 3(0 0 


I 2 \ I 22 I23 


c 2 (t) 


-c 2 (i) Ci(t) 0 


1 


. C 3(0 . 


where n is the SSF orbital rate of 0.0011 rad/ sec and the components of 
the nadir vector c(<) are given by the following expressions: 

ci (t) = -sin(9 2 (t))cos(9 3 (t)) 

c 2 (t) = cos(6i(t))sin(0 2 (t))sin(6 3 (i)) + sin(9 1 (t))cos(9 2 (t)) (3.4) 

c 3 (t) = —sin($i(t))sin(8 2 (t))sin(d 3 (t)) + cos(9i(t))cos(9 2 (t)) 

The attitude kinematics express the relation between the absolute an- 
gular velocities and the body-axis angular rates with respect to LVLH, as 
follows: 

$(t) = f (0(<)M<) + q (3.5) 


where 


f(0(O) = 


1 

cos(9 3 (t )) 


1 -cos(9i(t))sin(9 3 (t)) sin(9i(t))sin(6 3 (t)) 

0 cos (<9i(t)) -sin(9i(t)) 

0 sin{9 1 (t))cos(9 3 (t)) cos{9 1 (t))cos{9 3 (t)) _ 

(3.6) 



0 



q = 

n 

, m = 

9 2 {t) 


0 


. W) . 


The CMG dynamics are expressed in terms of the body-axis components 
of the CMG momentum and the external control torque u(<) applied by the 
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ACMM system as follows: 

0 — U>3 (t) U>2 (t) 

h(<) = u(t) - u> 3 (t) 0 -wi(<) h(t) (3.8) 

The general form of the external disturbance torque used in this study 
is as proposed in references [1,3]. That is: 

Wi(t ) = (Bias)i + Aisin(nt) 4- J3,sirc(2n<) (3-9) 

where i — 1, 2, 3 for roll, pitch and yaw, respectively. 



Chapter 4 

Derivation of the Linearized 
Equations of Motion 


The equations of motion presented in the previous section are highly non- 
linear and therefore not suitable for control system design. In order to allow 
use of linear control system design techniques some form of linearization 
must be performed. There have been successful applications of the small 
angle linearized and the roll/yaw and pitch decoupled equations of motion 
for the attitude control and momentum management of the SSF [1,4]. How- 
ever, for configurations resulting to large angle TEAs the applicability of 
this approximation will be extremely limited [3]. An alternate approach in 
resolving this problem would be the linearization of the equations of motion 
around an arbitrary TEA, without decoupling the pitch from the roll/yaw 
equations of motion. This is justified on the basis that the large angle TEA 
configurations of the SSF result from the misalignment of the control axes 
with principal body axes. Therefore, neglecting the cross-inertia terms for 
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large angle TEA configurations would result in highly erroneous linearized 
models. 

4.1 Determination of the Torque Equilibrium 
Attitude 

In order to perform linearization of equations (3.1), (3.5) and (3.8) around 
an arbitrary TEA, it is necessary to estimate the average TEA of a given 
configuration over an orbit, after equilibrium conditions have been reached. 
Assuming that (0 \ , 0 a 2 , 6' z ) represent the TEA (or the average value of TEA) 
at steady-state for a given configuration, these angles can be calculated 
using the following set of nonlinear equations: 

r gy ro (w*) + r gg {0 *) + Bias = 0 (4.1) 

where u>* is the absolute angular velocity of the spacecraft evaluated at 8 " , 
and 

f(0*) w* + q = 0 (4.2) 

where the function f(0) and q are defined in the previous section. The 
latter of these two system of equations can be used to solve for u>* in terms 
of O' obtaining the following relations: 

\ 

—nsin{ 8 3 ) 

<jj' = — ncos(0j)co.s(03) (4-3) 

nsin(6{)cos(6l) 

Substitution of these expressions into equation (4.1), reveals a system 
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of three nonlinear equations which, for a given inertia matrix I, could be 
solved to estimate the average value of the TEA. 


4.2 Linearization of the Attitude Kinemat- 
ics 

The attitude kinematics, expressed by equation (3.5), can be linearized 
around (0 * , w*) obtaining the following linearized attitude kinematics: 


0(t) = 


cos(O^) 


t 

0 


n 

< 

0 

S9i(t) + 

—rising) 


_ -n(cos(^)) 2 _ 


0 


$ 03 (t) + 


cos(^) 


— cos(^i)sin(^3) 

0 

Sl)\ (<) + 

cos(0J) 

0 


sin(9l)cos(9l) 


Su 2 (t) +(4.4) 


sin($l)sin( 6 3 ) 

— sin(dl) 
cos{9\ )cos(# 3 ) J 

where, for i = 1 , 2, 3, 


Su 3 (t) 


66i(t) = 6i(t) - 0 : 

SiOi(t) = u>i(t) - UJ* 


(4.5) 

(4.6) 


represent the perturbations of the attitude and of the absolute angular 
velocities from the TEA and from the angular velocities corresponding to 
the TEA, respectively. 
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4.3 Linearization of the Gyroscopic Torque 

The gyroscopic torque term, given by equation (3.2), can be linearized as 
follows: 

M) = 

where each of the partial derivatives in the above equation, evaluated at 
0(f) = 0*, depends upon the elements of the inertia matrix, I, and the 

TEA vector, 9 m . 

The partial derivatives of the gyroscopic torque term, with respect to 
the absolute angular velocity of the spacecraft, axe given by the following 
expressions: 

—Ii2sin(6i)cos(8z) — Ii3Cos(9l)cos(9 3 ) 

2Ii3sin(9 j) + h3Cos(8^)cos(9 3 ) + (In — l33)sin(9i)cos($%) 

h3sin(9‘)cos(9 3 ) — 2Ii2sin(9 3 ) -f (In — hi)cos(8\)cos(8 3) 

(4.8) 

— 2/23 COS (#1) COS (#3) — I\3sin(8 3) + (I33 — hi)^ n (^\) cos (^3) 

Ii2sin(0i)cos(OZ) + h3sin(9 * 3 ) 

2I\2COs(8\)cos(8 3 ) — J 13 .sm( 0 i)co.s( 03 ) + (In — hi )^ n (^) . 

(4.9) 


du>3 


I 12 sin(9 3 ) - 2I 2 3sin(8* 1 )cos(9;) + (I 22 ~ h3)cos(8l)cos(9 3 ) 
- I 12 cos(9l)cos(9 m 3 ) + 2Ii3sin(9l)cos($l) + (I 33 - In)sin(8' 3 ) 
Ii 3 cos(9\)cos(8 3 ) — h3sin(9 3 ) 


dr. 


9V T 


duJ2 


—n 


dr. 


m = -n 


dui 


(4.10) 
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4.4 Linearization of the Gravity Gradient 
Torque 

The gravity gradient torque term, given by equation (3.3), can be linearized 
as follows: 

A r„(0 = + %tm + §£* («)} < 411 > 

where each of the partial derivatives in the above equation, evaluated at 
0(t} — 8* , depends upon the elements of the inertia matrix, I, and the TEA 
vector, 8 *. The elements of the vectors containing the partial derivatives of 
the gravity gradient torque with respect to the spacecraft attitude angles 
axe given by the following expressions: 

(— I12A1 + I13B1 )Cl + {I22 ~ ^ 33 )(-^i — -®l) — 4I23A1B1 

(In — ^33)^1 Ci + I\2(B\ — Aj) + 2/13 Ai B\ 

{(In - In)Bi + I23 Ax) Ci + In(A\ - B\) + 2 / 12 A 1 B 1 _ 

( 4 . 12 ) 


dr. 


99 


dh 


= —3 TV 


— il 2 (A 2 Ci + C2B1) + I\3(B2C\ — C2A1) 
Iu(BiC 2 - C1A2) - l33(A 2 Cl - B1C2) 
Iu(AiC 2 + B2C1 ) — I22(AiC2 + C1B2) 

+-f 22 (Al- 4. 2 — B1B2 ) 

+/ 12 (BiB 2 — AiA 2 ) 

+/23(-Bl^2 — CiA 2 ) 


dr. 


99 


de 2 


= — 3 n 
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® T 99 

ee 3 


where 


—I 33 {B\B2 -f A\A2) — 2I 23 (A 2 Bi + A 1 B 2 ) 
+2I 13 (C 1 C 2 +A 2 B x ) + I 23 {A*C 2 + C 1 B 2 ) 
-\-I\ 3 (A\A 2 — B1B2) + 2 Ii 2 {A\B 2 — C1C2) 


(4.13) 


-3 n 2 


(I 22 — h 3 ){A\A 3 — B\B 3 ) — Iu{A 3 Ci — B\C 3 ) 
(I 33 — In)(A 3 C\ — B\C 3 ) — I\ 2 {A\A 3 — B\B 3 ) 
(In — 122 X ^1 C 3 + B 3 C\) -f I 23 {B\C 3 — A 3 C \ ) 

—2I 23 (A 3 B\ + AiB 3 ) — I\ 3 (B 3 C\ + A\C 3 ) 
+2I\ 3 (A 3 B\ + C\C 3 ) + I 23 (A\C 3 + B 3 C\) 
+/ 13 (B 1j B 3 + AiM) + 2Iu(A 1 B 3 - C,C 3 ) 


(4.14) 


Ai = co5(^)5m(^2)-sm(^3) + sm(^i) co5 (^2) (4-15) 

B\ = sm(0J)sm(02) s * n (03) ~ cos(0j )cos(#2) (4-16) 

Ci = sin(^)cos(^) (4-17) 

A 2 = 5in(^ )cos(#2)^ n ($3) + co.s(0i) sz ' n (^2) (4-18) 

B 2 = -cos(0I)cos(02)-sm(0 3 ) + sm(0J)sm(02) (4-19) 

C 2 = C05(^2) co<s (^3) (4.20) 


yl 3 = sin(0l)sin(6 2 )cos(9l) (4-21) 

B 3 = -cos(^i)6in(02)co5(03) (4.22) 

C 3 = -sm(02)s*’ n (03) (4-23) 
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4.5 Linearized System Equations 

The above equations can now be assembled to form the overall linearized 
dynamics for the rigid-body motion of the SSF in the vicinity of an arbitrary 
TEA. 

4.5.1 Orbital Dynamics/Kinematics Equations 

The following three dimensional vector definitions are introduced for the 
partial derivative terms of the external torques to simplify the mathematical 
form of the final equations: 

T* 

syr 


Each of the components of the above vectors correspond to the roll, pitch 
and yaw components of the external torques. Following the standard formu- 
lation [1,3, 4, 6], the following state vector is defined for the orbital dynamics 
and the kinematics, to be used for control system design purposes: 

xi(<) = [&*(*), #»(<). M*), ^ (0. 6W*)) T ( 4 - 26 ) 

In order to eliminate <$uj 2 (t) and substitute 9 2 (t ) in the system equations, 
the kinematic relation for 9 2 (t) is substituted in the orbital dynamics. The 
overall linearized orbital dynamics /kinematics then take the following form: 

Ei Xi(<) = Ai Xi(<) + Bi u(<) -|- B 2 w (t) (4-27) 

where Xi(t) is given by equation (4.26) and the matrices Ej, Ai, Bi and 
B 2 are as follows: 


dr. 


gyr 


doji 


dl r. 


99 


dOi 


(4.24) 

(4.25) 
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I n cos{9\) InCOs($Z) Ii3 «*(«;) + J 12 «n(*I) 
l2icos{0i) I 2 2COs{9^) I23 COs(6^) + l22sin(9i) 

Iz\COs{6\) I 32 co ^{$ 3 ) IssCOS^Bi) l32sin(0i) 

00 0 

00 0 

00 0 

1 0 0 

0 1 0 

0 0 1 

0 0 /i 2 nsm(#3) 

0 0 l22nsin{9^ ) 

0 0 h 2 nsin{9l) 

t jV(i)«mW) iOW?) 

r„V(2)co S (0;) r^(2)cos(et) 
r„V(3 )«»W) r„V(3 )««(«;) 

COs(0j) — 5m(^3)c05(^J) 

0 cos(9f) 

0 sm(9j)cos(9^) 





\ 
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Tgyr(l)cOs(6l) + Tg yr (l)sin(Ol) 
T gyr(2)cOs(61) + T* yr (2)sin(d;) 

r z gyT { 3WflJ) + i(3Wfl|) 
0 
0 
1 


Tgg(l)c°s(6i) 

Tg g (2)c°s(e;) 

Tg 9 (3)cos(0i) 

0 

0 


— ncos{8\ )cos(8^) 

Tgg(l)cOs(0;) 

T* g (l)cos(e*) + Tg g (l)nsin(8z) 

T gg( 2 ) COS (0l) 

Tgg{2) c os($l) + Tg g (2)nsin(8l) 

Tgg(Z)C0S(6 J) 

T gg(3)cos(0i) + T gg (3)nsin(8l) 

0 



0 0 
0 n6m(0j)5m(03) 


1 0 0 0 0 0 
0 1 0 0 0 0 
0 0 1 0 0 0 
0 0 0 0 0 0 
0 0 0 0 0 0 
0 0 0 0 0 0 


— Bi 


( 4 . 29 ) 


, B 2 = 


( 4 . 30 ) 
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4.5.2 CMG Dynamic Equations 

In addition to the above linearizations, the CMG dynamics can also be 
linearized as follows: 

0 — W3 U?2 

h (t) = u(<) wj 0 -wj h(t) (4.31) 

— IDj 0 

where the vector h(<) represents deviations of the CMG momentum from 
equilibrium conditions and where the vector u>* is given by equation (4.3) 
in terms of the TEA. As both the CMG momentum and its integral will be 
used in the feedback controller, the following state vector is defined for the 
CMGs: 

x 2 (i) = [hi (0,M0,M0,/ MO,/ MO,/ M0J T (4-32) 

Then the standard state-space form for the CMG dynamics becomes: 

x 2 (0 = A 2 x 2 (0 + B 2 u(i) (4.33) 

where x 2 (t) is given by equation (4.32) and the matrix A 2 is as follows: 

0 nsin($i)cos(6^) ncos(8i)cos(9 3 ) 

-nsin(B m 1 )cos(9l) 0 -na:n(0J) 

ncos(6l)cos(0l) nsin(B^) 0 

1 0 0 

010 
0 0 1 



CHAPTER 4. LINEARIZED EQUATIONS 4-11 

0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 

4.5.3 Disturbance Rejection Filter Equations 

As the frequencies of the external disturbances are known, the following 
disturbance rejection filters (DRFs) can be used for asymptotic rejection 
[ 1 ]: 


di(t) -f- n 2 ot\(t) — /ii(f) 

(4.35) 

Q 2 (f) + n 2 a 2 (t) = 0 2 (t) 

(4.36) 

a 3 (t) + n 2 a 3 (t) = d 3 (t) 

(4.37) 

ft(f) + (2n) 2 /M*) = M0 

(4.38) 

$i(f) + (2n) 2 /? 2 (f) = 9 2 (t) 

(4.39) 

fo(t) + (2 n) 2 fo(t) = 9a(t) 

(4.40) 


where, ctj(f), are the filter states with arbitrary initial conditions. It 

is worth mentioning that the above filters are for attitude hold, only. If 
momentum hold is desired then a separate set of DRFs should be designed 
with inputs hi(t), h 2 (t ), 0 3 (f). The reasons for not using 6i(t) or h 3 (t ) in 
the DRFs are given in reference [1]. The above equations can be expressed 
in the following state-space form: 

x 3 (i) = A 3 x 3 (t) + A 4 xi(t) + As x 2 (f) 



(4.41) 
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where X3(<) is given by: 


X3(<) = [di(i), C*l(0? ^(t), a 2(0» a 3(t)i 

Pi ( t ), Pi (t), Pl(t), Pi{t), Pz(t), &(<), f 

and the matrices A 3 , A 4 and A 5 are as follows: 

0 -n 2 0 0 0 0 

1 0 0 0 0 0 

0 0 0 -n 2 0 0 

0 0 1 0 0 0 

0 0 0 0 0 -n 2 

0 0 0 0 1 0 

Aa_ 000000 
000000 
000000 
000000 
000000 
000000 


(4.42) 
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0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 -4 n 2 0 0 0 0 

1 0 0 0 0 0 

0 0 0 —4 n 2 0 0 

0 0 1 0 0 0 

0 0 0 0 0 —4n 2 

0 0 0 0 1 0 


a 4 


000000000000 

000000000000 

000000000000 

000000000000 

001000000100 

000010000001 


100000100000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 


( 4 . 43 ) 


( 4 . 44 ) 


( 4 . 45 ) 
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The system of equations (4.27), (4.33) and (4.41) can now be trans- 
formed into the following state-space form: 

x(t) = Ax(<) + Bu(<) + Gw(<) (4.46) 

where the 24th order state vector is: 

x i(0 

x(<) = x 2 (f) ( 4 - 47 ) 

. X 3(*) . 

and the system matrices A (24 x 24), B (24 x 3), and G (24 x 3) are: 


E^Ai 

0 

0 


0 

a 2 

0 

(4.48) 

A 4 

As 

A 3 


[ E r 

*Bi 

1 



B = B 2 (4-49) 

0 



(4.50) 


with the rest of the submatrices previously defined and'dependent upon 
the specific value of the estimated TEA vector. By setting 8* = 0 for 

i = 1,2,3 and Lj = 0 for i ± j, in equation (4.46), the small-angle 
linearized and the roll/yaw and pitch decoupled SSF models are obtained, 
respectively. 


Chapter 5 

Design of the Attitude 
Control/Momentum 
Management System 


The dynamic system described by equation (4.46) is used to design an 
ACMM system. The disturbance rejection filters proposed in reference 
[1] are used to accomplish either attitude or momentum hold, depending 
upon the desired mode of operation. As mentioned earlier, however, only 
attitude hold is examined in this study. The results presented in this paper 
axe equally applicable to the momentum hold mode, as the only change 
occurs in the input variables to the DRFs. 
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5.1 Control System Architecture 

The overall architecture of the investigated ACMM system is shown in Fig- 
ure 5.1. The controller is composed of two gain- scheduled control laws. One 
for attitude control and momentum management, and the other for rejec- 
tion of the environmental disturbances. Both control laws are scheduled 
using the estimated average TEA as the scheduling variable. The input 
signals to the ACMM controller are the angles and the angular rates, the 
CMG momentum and the momentum build-up, i.e. the momentum inte- 
gral, and an estimate of the average SSF TEA based upon estimates of 
the available mass properties. The input to the DRFs are the angle and 
momentum vectors. However, as at most three independent variables can 
be simulaneously controlled, either the vector (hi(0> 0 2 (t), 0 3 (t)) or the vec- 
tor (M<), MO, 03(0) is used M in P ut to the DRFs ’ de P endin s u P° n the 
desired operating mode of the ACMM system. 

The use of the TEA estimate in the control law is necessary because of 
the inadequacy of the linear control system design tools to account for the 
controller gain variations with changing inertia matrix elements. The static 
calculation involved in obtaining the TEA estimate could be performed off- 
line, with table look-ups used during actual operation. It is assumed that 
a separate module provides on-line or tabulated estimates of the SSF mass 
properties. A saturation element is included as part of the simulated closed- 
loop system, because of the possible differences in the commanded and 
actual control signal that could be caused by high demanded control torques 
(the control torque saturation is set at 150 lb f , whereas CMG saturation is 
set at 20,000 ft - Ibj - sec). 
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The attitude/momentum control law used in this study is as follows: 

u (t) = -K(0*)x(t) (5.1) 

where the state vector x(t) is given by equation (4.47) and the 3 x 24 gain 
matrices, K, are determined for various TEAs, using the LQR algorithm 
described in the following section. 

5.2 LQR Design with Pole Placement 

Because the full state considered in equation (4.47) is available for mea- 
surement and because a time-optimal control law of minimal complexity is 
desired, the LQR is appropriate for use in obtaining the constant gains of 
the ACMM control system. The various linear designs are performed using 
the dynamic system described by equation (4.46), with the elements of the 
system matrices evaluated at specific values of the TEA, corresponding to 
various assumed configurations of the SSF during the assembly sequence. 

Assuming the state-space form of equation (4.46) with an n-dimensional 
state vector, an m-dimensional input vector and without the external dis- 
turbance term, the LQR problem can be stated as the minimization of the 
following functional: 

J = I (x r (f)Qx(t) + u T (t)Ru(t))dt (5.2) 

Jo 

where the weight matrices Q and R are (n x n)-dimensional non-negative 
definite and (rn x 7 n)-dimensional positive definite symmetric matrices, 
respectively. The full-state feedback control law that minimizes the above 
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functional is given by: 

u(f) = — Kx(f) = — R _1 B r Px(f) (5.3) 

where K is the feedback gain matrix and P, a (n x n)-dimensional non- 
negative definite symmetric matrix, is the solution of the following Riccati 
equation: 

PBR _1 B r P - PA - A t P - Q = 0 (5.4) 

where the pair (Q, A) is assumed detectable. 

It is usual practice to arbitrarily choose the matrices Q and R, resulting 
in a certain gain matrix K and the associated stable closed-loop poles. 
However, these closed-loop pole locations, though in the left-half plane, 
may or may not be desirable. The problem of selecting the matrices Q 
and R such that the resulting closed-loop pole locations fall in a desired 
sector, has been widely investigated. This is the so-called LQR problem 
with pole placement and it has most recently been addressed by Sunkel et 
al.[6]. Throughout this study the LQR with pole placement suggested in 
the above reference has been used with some modifications in the solution 
approach of the Riccati equation. The algorithm proposed by Sunkel and 
Shieh [6], allows the choice of the weight matrix R and the desired sector 
within which the closed-loop poles are to be located. The feedback gain 
matrix K and the weight matrix Q which place the closed-loop poles in the 
chosen sector are then determined. 

One of the numerical issues encountered in using the LQR problem for 
the SSF with large TEAs is the solution of the associated Riccati equation. 
As the order of the Riccati equation to be solved in this study is fairly high 
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- the closed-loop system is 24th order, however, the Hamiltonian formed 
is 48th order - use of the algorithms existing in the various control system 
design packages resulted in very inaccurate solutions. Therefore, the Riccati 
equation involved in the LQR/Pole placement algorithm was solved using 
the matrix sign function [6], with the following expression used for the 
evaluation of the matrix sign function: 

sign(A) = hm S (k) (5-5) 

where S(fc) is evaluated using the following recursive expression: 

S(Jfc + l) = S(fc)[5I + 10S 2 O) + S 4 (fc)][I + 10S(*) + 5S 4 (fc)] -1 (5-6) 
with S(0) = A. 

The LQR based control gains have been obtained for various assumed 
SSF configurations. The inertias for two of the configurations analyzed are 
given in Table 1, along with the corresponding estimated values of the TEA. 
The Assembly Flight #3 (MB3) corresponds to a highly non-symmetric SSF 
configuration, resulting in a large value of pitch TEA. In contrast the AC 
configuration corresponds to a highly symmetric configuration with small 
TEA. The numerical values of the gains for these two configurations can 
be obtained from the figures presented in the following pages. For all of 
the linear designs the closed-loop pole locations are chosen to be within 
the sector shown in Figure 5.2, as it has been the experience of the authors 
that containing the closed-loop poles in this region reveals relatively ’’good” 
system performance. 
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Figure 5.2: Pole Assignment Sector for the LQ Pole Placement Algorithm 
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Table 1. Space Station Freedom Configuration Parameters 


Parameters 

MB3 

AC 

Inertia ( slug — ft 2 ) 



I u 

23.22 x 10 6 

50.28 x 10 6 

1 22 

1.30 x 10 6 

10.80 x 10 6 

1 33 

23.23 x 10 6 

58.57 x 10 6 

^12 

O 

i-H 

X 

CO 

CM 

O 

O 

! 

-0.39 x 10 6 

I 13 

0.477 x 10 6 

0.16 x 10 6 

123 

-0.011 x 10 6 

0.16 x 10 6 

TEA ( degrees ) 




1.3 

0.19 

02 

-44.7 

-5.0 

oi 

2.1 

0.56 


5.3 Gain-Scheduling of the Linear Gains 

Gain-scheduling allows nonlinear control system implementation of linear 
control laws, which have been designed using any linear design method. 
A gain-scheduled control law is an open-loop adaptive control system [7]. 
Even though feedback control is used, the adaptation is open-loop and 
the system performance may suffer. A continuous gain-scheduled control 
system is designed by simply fitting the various LQR gain matrices to a 
certain nonlinear function. 

In view of the linearized system model of equation (4.46), the linear 
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control law obtained by repeated application of the LQR algorithm is: 

u (t) = K,x(f) (5.7) 

where i = 1,2,.. is the number of SSF configurations for which the LQR 
algorithm has been applied. Second order polynomials were used to sched- 
ule the elements of the gain matrices K^. The behavior of the fitted gains 
as a function of the TEA are depicted in Figures 5.3 through 5.25. 

As a result of the nonlinear nature of the gain-scheduled controller and 
of the open-loop system, there is no systematic procedure to determine the 
performance and the stability characteristics of the nominal SSF config- 
uration, other than to exhaustively simulate the system for a variety of 
scenarios and for sufficiently large number of orbits. Additionally, the ro- 
bustness properties of the controller can only be investigated via numerical 
simulations. The results of these studies are presented in the following 
section. 
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Figure 5.3: K(i,l) Gain Variation with the TEA. 
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Figure 5.7: K(i,5) Gain Variation with the TEA. 
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Figure 5.10: K(i,8) Gain Variation with the TEA. 
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Figure 5.11: K(i,9) Gain Variation with the TEA. 
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Figure 5.12: K(i,10) Gain Variation with the TEA. 
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Figure 5.13: K(i,ll) Gain Variation with the TEA. 
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Figure 5.14: K(i,12) Gain Variation with the TEA. 
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Figure 5.16: K(i,14) Gain Variation with the TEA. 
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Figure 5.17: K(i,15) Gain Variation with the TEA. 
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Figure 5.18: K(i,16) Gain Variation with the TEA. 
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Figure 5.19: K(i,17) Gain Variation with the TEA. 
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Figure 5.20: K(i,18) Gain Variation with the TEA. 












Figure 5.21: K(i,19) Gain Variation with the TEA. 
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Figure 5.22: K(i,20) Gain Variation with the TEA 
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Figure 5.23: K(i,21) Gain Variation with the TEA. 


PITCH TEA (DEGREES) 




5-31 



Figure 5.24: K(i,22) Gain Variation with the TEA. 









SPACE STATION RIGID BODY DYNAMICS LQR/POLE PLACEMENT GAINS 

INERTIA MATRIX VARIATION FROM MB 3 TO AC CONFIGURATION 


CHAPTER 5. DESIGN OF THE AC /MM SYSTEM 


5-32 



(ava/ai-u) samiNDvw nivo 


Figure 5.25: K(i,23) Gain Variation with the TEA. 
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Figure 5.26: K(i,24) Gain Variation with the TEA. 





Chapter 6 

Simulation Results 


The primary purpose of this section is to present some numerical simula- 
tions to substantiate the following two claims: 

(1) Certain SSF configurations characterized by large cross-inertia terms, 
such as the earlier flights in the assembly sequence or during some of 
the MRMS operations [2], require special attention when designing 
linearized control systems. 

(2) Because of the slowly time-varying nature of the SSF dynamics, a 
gain-scheduled control system provides "good” closed-loop perfor- 
mance, even for fairly large variations in the SSF inertia matrix. 
This results in less complex control laws as compared to alternate 
adaptive schemes requiring on-line identification, though, the system 
performance may suffer if the uncertainty in the inertia matrix is very 
high. 
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Figure 6.1: Aerodynamic Torque Profiles. 
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Figure 6.2: Roll Axis Angle Response for MB3 Configuration 




SPACE STATION CONTROLLED RIGID BODY DYNAMICS (NONLINEAR SIMULATOR) 


6 -' 



100 200 

TIME (MINUTES) 





SPACE STATION CONTROLLED RIGID BODY DYNAMICS (NONLINEAR SIMULATOR) 


CHAPTER 6. SIMULATION RESULTS 


6-5 



Figure 6.4: Roll Axis Control Torque Response for MB3 Configuration 
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Figure 6.5: Roll Axis Angle Momentum Response for MB 3 Configuration. 
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Figure 6.6: Pitch Axis Angle Response for MB3 Configuration. 
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Figure 6.7: Pitch Axis Angular Rate Response for MB3 Configuration 
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Figure 6.8: Pitch Axis Control Torque Response for MB3 Configuration. 
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Figure 6.11: Yaw Axis Angular Rate Response for MB3 Configuration. 
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Figure 6.12: Yaw Axis Control Torque Response for MB3 Configuration. 
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Figure 6.13: Yaw Axis Angle Momentum Response for MB3 Configuration. 
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All of the simulations subsequently presented in this section have been 
performed using the nonlinear simulation model for the SSF presented in 
section II of this paper. To demonstrate the first claim, two sets of simula- 
tions are performed using the MB3 configuration. In both cases the LQR 
synthesis algorithm with pole placement is used and the closed-loop poles 
are placed within the boundaries drawn by the two 45 deg lines and the 
— 0.5n vertical line, as shown in Figure 5.2. For the first set of simulations 
it was assumed that even though an accurate estimate of the TEA will 
be available, i.e. the TEA estimate will be obtained using the full inertia 
matrix, the controller design will be based on a model that ignores the 
cross-inertia terms. Therefore, these simulations utilized a linear control 
system designed using the model described by equation (4.46), with the 
cross-inertia terms set to zero. However, the actual values of the TEA, as 
given in Table 1, were used in the evaluation of the various non-zero matrix 
elements of the system dynamics, i.e. the LVLH approximation was not 
made. The second set of simulations were performed using a linear control 
system based on the fully coupled roll/pitch/yaw motion as described by 
equation (4.46), and with the same values of the TEA as before. In both 
sets of simulations, the same initial conditions and environmental distur- 
bance terms were used. 

The simulation results for these two cases are shown in Figures 6.1 
through 6.13. These studies clearly demonstrate the error introduced in 
arbitrarily discarding the cross-inertia terms when the analyzed configura- 
tions correspond to large angle TEA. The transient response of the closed- 
loop system becomes unstable within less than one-tenth of an orbit, unless 
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the cross-inertia terms of the linearized system model used in the control 
synthesis are included. The controller response with the cross-inertia terms 
included, is not perhaps the best possible, as further fine tunning of the 
desired closed-loop pole locations could possible reveal transient responses 
with less excursions. 

Demonstration of the second claim proved to be more time consuming 
and computationally more intensive. In order to arrive at a gain-scheduled 
controller, several linear designs were performed for various SSF configu- 
rations with TEA ranging form —45 deg to LVLH. These configurations 
correspond to inertia matrix variations of up to 1600%, for some of the 
cross-inertia elements. The gain-scheduled closed-loop dynamics were sim- 
ulated for 130 orbits while the inertia matrix was allowed to change in a 
stair-case fashion from the MB3 configuration to the AC configuration, as 
shown in Figures 6.14 through 6.19. The results of these simulations are 
presented in Figures 6.20 through 6.30. With the exception of the bias 
terms the external disturbances were kept the same, as the inertia varia- 
tions were made. These transient responses indicate that a gain- scheduled 
controller, whether it is in table look-up form or a continuous polynomial 
as presented in this study, is sufficient in handling large changes in the SSF 
inertia matrix, including configurations that are characterized by large an- 
gle TEA. This is primarily attributed to the slowly varying dynamics of the 
SSF, with characteristic time constants on the order of tens of minutes. 
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Figure 6.14: 1(1,1) Inertia Variation from MB3 to AC Configuration. 
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Figure 6.15: 1(2,2) Inertia Variation from MB3 to AC Configuration. 
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Figure 6.17: 1(1,2) Inertia Variation from MB3 to AC Configuration. 
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Figure 6.18: 1(1,3) Inertia Variation from MB3 to AC Configuration. 
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Figure 6.19: 1(2,3) Inertia Variation from MB3 to AC Configuration. 
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Figure 6.20: Roll Axis Angle Response for Inertia Variations from MB3 to 
AC Configuration. 
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Figure 6.21: Roll Axis Angular Rate Response for Inertia Variations from 
MB3 to AC Configuration. 
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Figure 6.22: Roll Axis Control Torque Response for Inertia Variations from 
MB3 to AC Configuration. 
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Figure 6.23: Roll Axis Momentum Response for Inertia Variations from 
MB3 to AC Configuration. 


TIME (ORBITS) 






CHAPTER 6. SIMULATION RESULTS 


6-27 



Figure 6.24: Pitch Axis Angle Response for Inertia Variations from MB3 
to AC Configuration. 
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Figure 6.25: Pitch Axis Angular Rate Response for Inertia Variations from 
MB3 to AC Configuration. 
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Figure 6.26: Pitch Axis Control Torque Response for Inertia Variations 
from MB3 to AC Configuration. 
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Figure 6.27: Pitch Axis Momentum Response for Inertia Variations from 
MB3 to AC Configuration. 
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Figure 6.28: Yaw Axis Angle Response for Inertia Variations from MB3 to 
AC Configuration. 
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Figure 6.29: Yaw Axis Control Torque Response for Inertia Variations from 
MB3 to AC Configuration. 
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Figure 6.30: Yaw Axis Momentum Response for Inertia Variations from 
MB3 to AC Configuration. 






Chapter 7 

Summary and Conclusions 


This study presents a systematic derivation of the SSF dynamics linearized 
around an arbitrary TEA. Linearized models which appeared in the lit- 
erature can be obtained from the presented model by eliminating certain 
variables or by setting them to specific numerical values. The usefulness of 
the fully coupled model in control system design for large angle TEA config- 
urations is demonstrated by designing linear control laws based on models 
with and without the contribution of the cross-inertia terms. The nonlin- 
ear simulations indicate that control laws for ACMM based on small angle 
linearized or pitch and roll/yaw motion decoupled models of the SSF may, 
under circumstances, result in unstable closed-loop systems, even though 
the linearized closed-loop systems appear stable and with good overall per- 
formance. Additionally, numerical robustness studies have revealed that 
the closed-loop system remains stable for up to ± 30% variations in the 
SSF inertia matrix, as previously reportes in the literature [1]. 

As the inertia properties of the SSF will be varying for a number of rea- 
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sons, such as the assembly sequence and MRMS operations during or after 
assembly complete configuration, it is desirable to allow smooth ACMM 
throughout such scenarios. From this study it is evident that a single 
set of gains will not allow the controller to stabilize the system. It has 
been therefore proposed to use a gain-scheduled controller with the TEA 
being the scheduling variable. The approach suffers, as estimates of the 
TEA must be available in order to update the gains, however, it compares 
favorably to alternate adaptive control techniques requiring on-line identi- 
fication. Some additional computer storage will be required to retain the 
second order polynomial coefficients in memory. Alternate implementa- 
tions, such as a tabular form of the gains, could be used to further reduce 
the storage requirements. Transient response simulations indicate that for 
inertia matrix element changes of up to 1600%, the gain-scheduled con- 
troller provides ’’good” ACMM. It is not the intention of this study to 
fine-tune the gain-scheduling for best performance and/or robustness. Ad- 
ditionally, no attempt has been made to determine the maximum variation 
in the inertia matrix that the gain-scheduled controller can tolerate. These 
and other practical implementation issues, such as controller complexity 
and storage requirements, remain to be explored. 
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Part III 


Controller Synthesis for 
Input-Output Tracking in 
Uncertain Nonlinear Systems 



Chapter 8 
Introduction 


One of the fundamental design requirements for feedback control systems 
is servo-tracking or the ability to regulate controlled outputs without any 
steady-state errors in the presence of disturbances and plant parameter 
variations. This problem has been extensively studied in the literature, 
especially in the context of linear-time-invariant (LTI) systems. It is worth 
emphasizing here that in the classic asymptotic servo- tracking problem, 
the one that has been treated extensively in the literature [e.g., Francis 
and Wonham (1975)], no particular attention is given to the initial track- 
ing error, though it is clear that the initial errors will be bounded as a 
consequence of closed-loop stability. In the design of such systems the in- 
ternal model principle proposed by Francis and Wonham (1975) plays an 
important role. The internal model principle states that the controlled 
output tracks a class of reference inputs with zero asymptotic error, if the 
reference is included in the stable closed-loop dynamics. For example, recall 
that no steady-state error occurs for step reference commands in a type 1, 
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stable, closed-loop system which has an integrator in the loop. Note that 
the integrator is the generator for step inputs. This type of a servo system, 
although is robust to plant parameter variations, is quite sensitive to vari- 
ations in the reference input. For example, in the case of a step input, if a 
small sinusoidal variation is added, the steady state error will no longer be 
zero. The reason is that the integrator is not the generator of the reference 
command anymore. Consequently, at steady-state the closed-loop system 
will exhibit a sinusoidal oscillation. In real-world problems, however, one 
has to deal with such situations, and the internal model principle alone can 
not resolve these issues. 

In order to deal with potential inaccuracies in reference command in- 
puts, a novel notion of tracking, called tracking in the sense of spheres was 
recently introduced by Jayasuriya et al (1984). The emphasis there is that 
servotracking should be viewed as the ability to follow a desired output 
within a specified error bound over its entire duration; that is, to obtain 
a closed-loop system that tracks reference commands and rejects a class 
of external disturbances, to the extent required by the error bound, in the 
presence of plant uncertainties. 

In asymptotic tracking, usually a desired output vector y 0 (^) € R l is 
specified and the actual output vector y (t) € R l is required to follow the 
desired output asymptotically. The idea of tracking in the sense of spheres, 
unl ike the asymptotic tracking case, requires that ||y — y 0 || be smaller than 
a prespecified real number. Here ||.|| denotes any Zrj, - norm though the 
specific norm we consider in the present study is the Z^o ■ norm. 
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Remark 1: (i) Let y,y 0 G L l p , then 

||y — yolU*, = m f x (J 0 |y«(0 - yoi(t)\ p dty /p , i = i,2,...l 

where p is a finite positive integer. 

(ii) If y,y<> G L l <x>i th en 

||y - yolU' = max[esssup \y,(t) - y«(<)|], i = l,2,...7. 

t t £R 

The notion of tracking in the sense of spheres is primarily based on 
the fact that a system with uncertainties cannot be forced to track a given 
output precisely although it may be possible to do so with a small tracking 
error. When the tracking specifications are posed in a quantitative manner, 
as is done in this study, it becomes necessary to deal directly with the 
uncertainties of the system. Unlike in the classic servotracking problem, the 
necessity of having to deal directly with the uncertainties in this precise 
tracking problem makes it difficult to obtain criteria and design algorithms 
that are both sufficient and relatively simple. 

In the literature the tracking problem in the sense of input-output 
spheres has been formulated by Jayasuriya et al (1984), Jayasuriya and 
Kee (1988), (1990) in terms of topological neighborhoods in normed func- 
tion spaces Too and L 2 . The main design criterion central to the method- 
ology of Jayasuriya et al (1984) can be stated as a quantitative pole place- 
ment (QPP) procedure for adjusting the size of a certain linear operator 
norm. In particular, the sphere specifications are satisfied provided a proper 
placement of eigenvalues can be found so that a certain closed-loop opera- 
tor characterized by these eigenvalues belongs to a set of linear operators 
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whose Loo or L 2 induced norms axe upper bounded by a real number de- 
termined by the radii of input-output spheres and the uncertainties in the 
system. However, it has been difficult to obtain such eigenvalue locations 
in a systematic way. 

In Jayasuriya et al (1984) a trial and error procedure was suggested for 
the Loo - problem. In a recent paper, Jayasuriya and Kee (1990) studied the 
problem further and gave additional insight to the QPP. In particular, they 
proposed a certain perturbation of the classic Butterworth pole patterns as 
a solution to the Loo , QPP problem. Jayasuriya and Kee (1988) also 
studied the QPP problem when the tracking specifications are in terms 
of L 2 - measures and the plant is of the Lure type. They gave a circle 
criterion for solving the latter QPP problem which was facilitated by the 
Hilbert space structure of the function space 

A formal procedure for selecting the eigenvalues for the QPP problem, 
especially when Loo - tracking measures are employed, is currently lacking. 
The present study is aimed at partially filling this gap. We say partially be- 
cause the class of nonlinear systems considered is limited to be of the Lure 
type. A Lure type system lends itself to an L 2 - formulation which can sub- 
sequently be adopted to an L^ - framework. Embedding the problem first 
in L 2 allows us to draw upon some classical results from nonlinear systems 
theory on one hand and to utilize H - minimization [Francis(1987)] for 

V < 

solving the QPP on the other hand. Loo - performance specifications are ex- 
tracted from the L 2 - formulation by resorting to an exponential weighting 
of L 2 - functions, yielding Loo - functions. 

Part III of this report is organized as follows. First, some preliminaries 
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axe given in chapter 9. The input-output tracking problem related to Lure 
type nonlinear systems is formulated in chapter 10. In chapter 11, an 
algorithm for minimizing the i/oo-norm is given, followed by a step by 
step design procedure for the nonlinear tracking problem in chapter 12. 
A simple second order system is used in chapter 13 to demonstrate the 
design methodology, followed by a two-degree-of-freedom (2 DOF) robotic 
manipulator in chapter 14. Summary and conclusions from part III of the 
report are presented in chapter 15. 



Chapter 9 

Some Preliminaries 


Some basic concepts from elementary functional analysis and notation are 
collected below. 

Lp [0 , oo ) space: is defined as the space of measurable functions with 
the norm defined as 

||f||„ = naax( f° |/,(<)| p dt) 1/p < oo, * = 1,2 

where f(<) € R k . ||.|| P denotes the L p - norm of the function (.) in the 
L k p - space. When it is understood that f (<) € R k we will simply use the 
notation f G L p [0,oo) instead of f G Lj[0,oo) throughout this paper. 

When p = 2 we have the Hilbert space L 2 specified by f G L 2 [0,oo) 
where 

||f|| 2 = mpc ( J \fi(t)\ 2 dty /2 < oo, t = 1,2, ...k 

When p = oo we have the Banach space Leo defined by f G Loo[0, oo), where 

|| f ||oo = max[ess sup |/j(t)|] < oo, i = l,2,...fc 
* te[o,oo) 
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L 2e [0, oo) denotes the extended L 2 [ 0, oo) space defined as 


ft 


with 


{/ : 9 — V| vfe», ||(/)t ||2 = {[ \m\ 2 dtf 2 < 00 } 


/f ” 1 


/ 

0 


0 < t < f 
t > T 


where T denotes a truncation (Desoer, 1975). 

Central to the present study is the novel concept of tracking based on 
topological notions first introduced by Jayasuriya et al (1984). This notion 
of tracking is referred to as tracking in the sense of spheres. Also used in 
the study are ideas drawn from the Hoo - optimal control literature and the 
multivariable circle criterion cited below. 


9.1 Tracking in the Sense of Spheres 

Conventionally, tracking is considered to be the following of a specified 
trajectory y(t) in an asymptotic sense. That is, the actual trajectory y(t) 
should approach the reference trajectory y a (t) as t — » 00 . As opposed to 
this classic definition of asymptotic tracking, the present work is based 
on a more stringent tracking requirement where the actual trajectory y (t) 
is forced to follow the reference y 0 (t) as closely as possible for all times 
t G [0, 00 ). This is the notion of tracking in the sense of spheres which is 
illustrated in Figure 9.1 and expressed in a mathematically precise fashion 
below. 

Definition 1 : A given output (input), y : T — > R p (r : T — > R m ) is 
said to belong to an output (input) sphere fi(y,y 0 ,A>) (ft(r, r o ,0 T i) ) of 
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radius > 0 (/? r « > 0) centered at y 0 : T ♦ R? (r© : * R ) if 

||y _ y 0 || < jl 0 (||r — r 0 || < /?„•), where ||.|| is any norm associated with the 
output (input) function space T = {y|y • T — ► F} (A = {r|r : T —> R}) 
and T is the time set. y 0 (r 0 ) is referred to as the nominal output (input). 

Definition 2: If for any input r in the input sphere ft(r, r c , A-,), the system 
output y stays in the output sphere fl(y,y 0 ,A), then, we say that the 
system tracks the reference y 0 in the sense of input-output-spheres of radii 
A. and fio, respectively. 


9.2 The Multivariable Circle Criterion 

The multivariable circle criterion (Cook, 1972) can be summarized as fol- 
lows: 

Let a nonlinear, time-invariant system have the equations 

x(f) = Ax(f) + Bu(<) (9.1) 

y(t) = Cx(<) 

where the nonlinear control u (t) = — Fy(f), F = diag[fi(t,y )] and otiy 2 (t ) < 
/i(*,v(0)v(0 < fty 2 (0» < A, *' = l,2,...,m. Also, x(*),u(0 and y(t) 

are real vectors, of dimension ti , m and tii respectively. 

Let a = diag(ai) and (3 = diag(fii) such that a < F < (3, and the 
m x m transfer function matrix H(s) be defined as 

H(s) = C(sl - A) -1 B = [Aj(s)]> i = l,2,...m, j = l,2,...m 
The resulting closed-loop system is described by the nonlinear, time varying 
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Figure 9.1: Illustration for Precision Tracking. 
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equation 

x(<) = { A - B (diag[fi(t, y(t))])C}x(t) . (9.2) 

Now define 

di(s) = j£i\hij(s)\) i = 1,2, ...m 

and 

<<(.) = 

and let 19 be a contour consisting of the imaginary axis from s = —jR 
to s = j R, together with a semi-circle of radius R in the right half of the 
complex s- plane. Defined below are the ith Gershgorin mean band and the 
critical disc. 

Definition 3: The band swept out by the discs with center hu($) and 
radius [d,(s) + d'(s)]/ 2 as the complex variable s traverses D by the disc 
with center h ti (s) and radius [d,(s) + dj(.s)]/2 is named the ith Gershgorin 
mean band. The disc with the two points (— l/a t ,0) and ( — 0) as 
diameter is called the ith critical disc, The ith Gershgorin band 

and the critical disc are illustrated in Figure 9.2. 

With the above definition, the multivariable circle criterion (Cook, 1972) 
can be stated as lemma 1. 

Lemma 1 ( The Multivariable Circle Criterion): Let (A,B,C) 

be a minimum-order realization of H(s) and let the ith )3ershgorin mean 
band have no point in common with the ith critical disc, and let the ith 
Gershgorin mean band encircle the ith critical disc N cl times clockwise as 
s goes once clockwise around D, these conditions remaining satisfied for all 
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Figure 9.2: Illustration for zth Gershgorin Band and the Critical Disc. 
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such contours D with radius R > R, and suppose 

V? =1 N c i = - P 0 ( 9 * 3 ) 

where P 0 is the number of poles of the system transfer matrix H($) lying 
in the closed right half s-plane. Then the system (9.2) is stable in the sense 
of Liapunov. 

Proof: See Cook (1972) 

Lemma 1 can be viewed as an extension of the Nyquist stability criterion 
to Lure type systems with the actual disc playing the role of the —1 point 
and the Gershgorin mean band playing the role of the Nyquist plot. 

The relationship between Liapunov stability and the L 2 stability is de- 
scribed in Desoer and Vidyasagar (1975) and can be stated as: If the 

nonlinear control u (t) in Eqn. (9.1) is of the form -Fy (t) + r (<), where 
r(t) is a reference input, then, the closed loop system is L 2 stable (i.e., 
y(t) £ L 2 whenever r (t) £ L 2 ) provided the conditions of lemma 1 are 

satisfied. 


Chapter 10 

Problem Formulation 


In this section, the - tracking problem for Lure type systems is for- 
mulated. The formulation consists of first showing the I<2e - stability of 
a certain exponentially weighted form of the actual system equations, to 
assure L 0 0 - bounded-input, bounded-output stability of the actual system, 
followed by an estimation of the tracking errors in terms of the - mea- 
sures. The main design criterion is shown to be a ’’minimization of the 
Hoo - norm of a certain linear closed-loop operator. The approach to the 
problem is based on fixed-point techniques and some results drawn from 
nonlinear systems theory, especially the multivariable circle criterion. 

10.1 Plant and Control Structure 

Consider the nonlinear, uncertain system in the state-space form 

x(<) = Ax(<) + Bu(f) + Gd(t) — f(y(t), 7 , t) (10.1) 

y = Cx(f) 


10-1 
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where x(<) € JT is the state, y(t ) E R? is the output, u (t) E R m is 
the control, d (t) E IF is the disturbance vector, t € T = [0, oo) and 
f(y(t), 7 ,t) : IF x IF x T IF is a nonlinear function. A, B, G, C are 
constant matrices of appropriate dimensions. 

The control objective is to guarantee Too - robust tracking in the sense 
of spheres. That is we require the output y (t) to be inside a sphere of radius 
00 with the nominal output y c (t) as center for every reference input r(t) 
lying inside a. sphere of radius /?„• with the nominal reference input r 0 (t) as 
center. In Jayasuriya and Kee (1988) a similar problem was studied but 
the input and output spheres were defined in terms of the Li - norm. 

The controller structure shown in Figure 10.1, is adopted for effecting 
tracking in the sense of spheres where S c denotes a servo-compensator, S 0 a 
nonlinear observer and S a a stabilizing compensator. State-space descrip- 
tions of each of the blocks S c , S 0 and S, are given below. 

10.1.1 Servo-Compensator S c : 


x c (t) = A c x c (<) + B c (r(t) - y(t)) (10.2) 


The plant dynamics and the servo- compensator dynamics lead to the aug- 
mented state equations 


X a (0 


A 0 
-B c C A c 


x(t) 

Xc(0 


B 

0 


u{t) 


0 G 


r(() 


f(y 

B c 0 


. d w . 


0 


-3 
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= A a x a (t) + B a u(<) + G a w(i) f a (y(i),^,f) 


y(0 = [c o] 


x(f) 

X c (t) 


C 0 x 0 (<) 


(10.3) 

(10.4) 


where r (t) £ R v is the reference input and (A C ,B C ) is a controllable pair 
with |AI — A c | = A° + Qa— i A° 1 + ... + QiA + ot 0 . In addition A c £ R pq pq ■, 
B c £ R? qXp and x c £ R pq . The purpose of including a servo- compensator 
in the control loop is to assure zero error at steady state when there is no 
input uncertainty. 


10.1.2 Nonlinear Observer S 0 : 

In order to estimate the plant states, we use the nonlinear observer S 0 
having the following state equation: 

i(() = Ax(<) + Bu(t) - t(y(t)i 0 + H(y(<) - y(f )) (10.5) 

y(t) = Cx(f) 

where x(f) £ R n is the estimated state and H is a constant observer gain 
matrix of order n x p. f 0 (y (t),<) is a nonlinear design function obtained 
by replacing the uncertain parameters 7 in f(y(t), 7 i 0 by their nominal 
values. Note that although the augmented state is ^x (t) x c (t)J , there is 
no need for an observer to estimate both x T (t) and since x^(<) is 

already available by the construction, given in Eqn. (10.2). 
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10.1.3 Stabilizing Control Law S s : 

The stabilizing control law S s consists of state feedback from x(f ) and x c (t) 
and a nonlinear control fio(y(f) ? f) P ven by 

u(t) = [Ki K c ] x a (t) + fio(y (*), *) ( 10 - 6 ) 

where Ki and K c are constant state feedback gain matrices. The nonlinear 
controller fio(y(0»0 is obtained by replacing the uncertain parameters 7 
in fi(y(<),7,<) by their nominal values. 

We make the following assumptions regarding the state-space represen- 
tation of Eqn. (10.1): 


(i) The triple (A,B,C) is controllable and observable. 

(ii) f(y(<),7,<) satisfies the matching condition : 

f(y(<)>7,<) = Bfi(y(t),7,f) 

fc ( y (*)»0 = Bfio(y,t). 

(iii) There exists a constant matrix N such that B = GN. 


Now, by collecting equations (10.1 - 10.6), we arrive at the following 
closed-loop equations 

z (t) = A,z(0 + B,w(t)-^i(y(<),7,0 
= A z z(t) + B,w(t) - B z (f>(y(t), 7 , t ) 
y(t) = C z z(t), 


(10.7) 
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where 



A BK C BKi 


0 

G 

Ill 

N 

◄ 

— B c C A c 0 

, B* = 

B c 0 


HC BK C A + BKi - HC 


0 

0 J 


fi(y(0>7><) - My (*)>*) 


c, = [C 0], ^a(y(0i7i0 = B 


0 

0 


^(y(0i7i0 




x(t) 

0 

(fi (y(t), 7, t) - fio(y (0, 0). 311(1 z ( t ) = 

Xc(t) 

N 

. x(0 . 


In the above closed-loop arrangement the design parameters are the 
matrices Ki,K c ,H and the nonlinear function fio* 

Since the control objective is Leo - tracking in the sense of spheres, 
we now proceed to establish BIBO stability of the exponentially weighted 
closed- loop system with an appropriate bound on Loo tracking error. In 
order to establish the results, we first prove L 2 - stability of the closed-loop 
system by exploiting the fact that Eqn. (10.7) is of the classic Lure form. 
The circle criterion gives conditions for L 2 _ stability. Once L 2 - stability 
is guaranteed, we conclude Loo stability of the closed-loop system by the 
exponential weighting technique. After proving Loo - stability, we provide 
estimates for the tracking errors in terms of the Loo * norm. One primary 
motivation for embedding the problem first in L 2 is to utilize - optimal 
control theory for obtaining the design matrices K = [Ki,K c ], and H. 
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10.2 X 2 - Stability 

Consider the closed-loop system shown in Figure 10.2, described by Eqn. 
(10.7). 

By defining H(s) = C z (sl — A Z ) - 1 B Z with h(t), the impulse response cor- 
responding to H(s), Eqn. (10.7) can be written in the following integral 
form. 

y(t) = J h(t — r)w(r )dr — jf h(t — r)</>(y(r ), 7 , t)c?t 

= v(t) / h(t - r)<^(y(r), 7 ,r)dr, ( 10 . 8 ) 

Jo 

where v(t) = / h(< — r)w(r)dr. 

Jo 

On multiplying both sides of Eqn. (10.8) by e at ,a G R 1 , we obtain 

e at y(t) = e at \(t) — f e a(t_T) h(< - r)e aT </>(y(r), 7 , r)dr (10.9) 

J 0 

which can be written as 

y B (<) = v„(<) - f h a(t - r)e aT </)(e" aT ya(r),7,T)dr , (10.10) 

v 0 

on substituting y a = e a ‘y( 0 >h a = e a *h(t), and v 0 = e at v(t) = f 0 h a (t — 
r)w n (r)dr with w n (t) = e at w(t). 

We assume that the nonlinearity </>(y(t), 7 , <) lies in the sector [a, /?], /3 > 
0 (Desoer, 1975). With a = diag(au), (3 = diag((3i ) and <p = diag(4> t ), i = 
l,m, where m is the dimension of <£, the sector condition becomes a.yf < 
yifc < (3 t y1 for i = l,m. We further assume that h a = e a< h(t) G Li flL 2 , for 
some a > 0. The L 2 - stability of the exponentially weighted system (10.10) 
can be obtained by the multivariable circle criterion and is summarized in 
the following lemma : 
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Lemma 2: If the Nyquist diagram uj H a (ju) = H(-a + jiv) satisfies 

conditions of the multivariable circle criterion, then, the system of Eqn. 
(10.10) is Z< 2 e stable (i.e., w n € L 2 e => Ya € Lie )• 

Proof: From Eqn. (10.10), we have 

y.(t) = v a (<) - jf h a (t - r)e aT ^(e" 0T ya(r),7, r)dr (10.11) 

which is equivalent to Eqn. (10.9) 

e ot y(t) = e at v(t) — J* e a(<-T) h(t — r)e aT (/>(y(r), 7 , r)dr 
= e ai I* h (t - r)w(r)dr - e a(t - T >h(* - r)e OT 

Jo Jq 

4>(y(T),i,T)dT (10.12) 

The block diagram corresponding to equation (10.12) is shown in Figure 
10.3. 

Since the output y a = e a< y(t) and «fty(<),7,0 G [«,£] with respect to 
y(t), it can be easily concluded that e at </>(y(<),7,t) = e at <f>(e a y a (<),7 >0 
belongs to the same sector [ot,(3] with respect to y a (t) because the latter 
simply involves a scaling by e at (see Figure 10.4). 

Since <f> a = e at 4> £ [a,/3], we can directly apply the circle criterion to 
determine the stability of system (10.10). Consequently when w n (<) £ L 2e 
and H a (ju) satisfy the conditions of the circle criterion or lemma 1, it 
follows that y a £ L 2e . 

Note that since H a (jw) = H(— a, + jw), the L 2e stability of the system 
given by Eqn. (10.10) is verified by replacing the standard polar plot U(ju) 
by H(— a + ju:) in the circle criterion. 
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Figure 10.3: Block Diagram of the Exponentially Weighted Closed-Loop 
System. 
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10.3 Zoo - Stability 

In what follows, we use the L 2 e * stability conditions of the exponentially 
weighted system (10.10) to establish the L qq stability of the actual closed- 
loop system represented in the integral form in Eqn. (10.8). 

Lemma 3: If the exponentially weighted system is L<it stable then the 

original system is L ^ stable. 

Proof: Recall that system (10.10) is 

y a(t) = Va(t) ~ £ h 0 (t-r)e“^(e' aT y a (r), 7 ,r)(ir (10.13) 

where the nonlinear element G [a,/3]. Since system (10.10) is L 2e stable, 
there exists some finite constant p such that 

||[y„(.)]r||2 < HIMOHI2 (10.14) 

where ||[y a (-)Wl2 denotes the L 2 norm of the truncated function [y a (-)]T- 
Using the fact that h a € L\ and 

\ a (t) = I h a (< — r )w n (r)dr (10.15) 

Jo 

and the Parseval Theorem [Desoer and Vidyasagar (1975)], we obtain 

ll(v.(.)M» < sup||H.0'«)III|[w»(-)]t||.. (10.16) 

Substitution of Eqn. (10.16) into Eqn. (10.14) yields 

|| [y a (-)]r 1 1 2 < />sup||H a 0'w)||||[w n (.)] T ||2 (10.17) 

u> 

and from the proof of the circle criterion [Desoer and Vidyasagar (1975)], 
we obtain 

llly-OMIr < >IIK(-)Wl 2 (10.18) 
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where 

i = \\\P-<x\l m = \\\P + Q W > 

H «(«) = (I + mH a (s))" 1 H a (s) 

70 = sup l|H t (icc?)l|. 

U) 

By comparing Eqns. (10.17) and (10.18), it now follows that 

P = (sup ||H.0 'k>)||)' 1 t-^— (10.19) 

substitution of which in Eqn. (10.14) yields 

||[y.(.)]r|| 2 < (sup||ff«aw)||r 1 — ^-||[v a (.)W| 2 . (10.20) 

uj 1 — t /0 

It can be easily seen that when w (t) £ L<x>, then v(<) £ Loo, since h(t) £ L\ . 
Moreover, 

IIMOMb 3 |||e"v(.)]r|| 2 = (jf l|e"'v(<)H ! *) ,/! 

r T r (p 2aT — i y/ 2 ' 

< sup||v(t)||(^ e 2o< d<) 1/2 = ||v(.)||oo (2a) 1 / 2 

< (2 a) _1/2 \/e 2o:r - l||h(.)||i||w||oo • (10.21) 

Since the nonlinear element e at 4>(e~ at y a (t), 7,<) belongs to the sector [a,/3] 
with respect to y a (t), we obtain 

||[e ot <5!)( e - ot ya(0,7,0]r!|2 < /3||[ya]x|| 2 (10-22) 

which becomes 

||[e“ , <£(e"“y.(f).7. ( )]rl|! < sup ll H .(j‘‘ J )ll)~ 1 , J° 7o 

((2 - l)||h(.)||, 


w||<410.23) 
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upon substitution of Eqns. (10.20) and (10.21). Taking the amplitude on 
both sides of Eqn. (10.8), we obtain 

|y(OI < |v(*)| + / \H t - T )\\ < P(y( T )^,T)\dT 

J 0 

< ||h||i||w||oo + e _at I* e o(t_T) |h(< — r)|e aT |<^(e _OT y a (T), 7, r)\dr 

JO 

< l|h||i||w|U + e— ||h(r)e"|| 2 ||[e"V(c-‘y.(<).-l'.<)H|2 

(from Schwartz Inequality) 

< IMiMloo + ||h(r)e OT ||2||[e a< </>(e- at ya(0,7,<)]T|| 2 (10-24) 

(since t > 0, a > 0). 

Now using Eqn. (10.23) in (10.24), the following is obtained. 

||y(OI|oo < [||h||i + ||h(r)e aT || 2 ^(sup|H a (ja;)|)" 1 

{(2a)" 1/ V e2aT - l}l|h(.)||i]||w||oo 

= Kk || w||oo , where Kk is a finite constant. (10.25) 

From Eqn. (10.25) it is now clear that if w (<) G Loo ( he., ||w(t)||oo < °°) 
then ||y(t)||oo < °o implying that the original system is BIBO stable. This 
completes the proof of Lemma 3. 

Lemmas 2 and 3 are now combined in theorem 1 below to give conditions 
for guaranteeing L 0 0 - stability of the closed-loop system (10.8). 

Theorem 1 : Given the nonlinear equation (10.8), with h(.) G L\ , e“'h(.) G 

Li nL 2 for some a > 0, and with the nonlinear element <f>( y(.), 7, .) belong- 
ing to the sector [o:,/ 3], (3 > 0, the closed-loop system (10.8) is Loo stable, 
provided that the Nyquist diagram u> i-> H a (ju;) = H(-a + ju>) satisfies 
the conditions of the multivariable circle criterion given in lemma 1. 

Proof: The proof of theorem 1 is immediate from lemmas 2 and 3. 
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10.4 Loo - Tracking 

From theorem 1 we note that the closed-loop stability in the L <*, - sense is 
established by considering the modified Nyquist locus H(-a+ju;) resulting 
from Eqn. (10.12). The transfer matrix H(-a 4- 5 ) can be thought of as 
representing a system equivalent to the modified state equations 

z a (t) = (al + A.r)z a (i) + B z w„(t) — B z <f> a (t) (10.26) 

y a(t) = CjZa(t) 

having the integral form 

y a (t) = C z e A " (1_T) B,w„(r)dr - 

C z / e An( ' _T) B z </) 0 (ya(r), 7 ,r)dr + C z e Ani z o (0) 

Jo 

where 

A n = al + A 2 (10-27) 

or the operator form 

y a (<) = C z rB 2 w n (t) - C z rB z <t> a (t)y a (t) + y a (0) (10.28) 

where y a (0) = C 2 e Ant z o (0) denotes the system initial conditions, and 

(rx)(<) = J e An(< - T )x(r)dr, 

<M0ya(0 - ^a(yo(0>T>0 £ l a i Pi- 

Remark 2: Since (A,B,C) is controllable and observable, the control 

gains K = [Ki, K c ] and H can be found so that A n (= (al-f A z )), a > 0, is 
Hurwitz. Moreover, under the assumption that (f>{y (t),7,t) in Eqn. (10.8) 


CHAPTER 10. PROBLEM FORMULATION 


10-16 


is sector bounded in [a,/3], and the circle criterion in lemma 1 is satisfied, 
we know that the system (10.28) is L 2e stable because also belongs to 

the same sector. Thus, when w n belongs to L 2e , the linear and the nonlinear 
operators T and <f> a (t), map L 2e back into itself. 

We now recall that our control objective is to effect L^ - tracking in the 
sense of spheres for the uncertain nonlinear system (10.8). So far we have 
established £«> stability of (10.8) and L 2e - stability of the exponentially 
weighted system (10.10). Because one of our goals is to specify the design 
matrices Ki,K c ,H by an - formulation, we first solve the tracking 
problem in L 2e to obtain H <*, - conditions on the linear closed loop operator 
(C*rB 2 ) and subsequently enforce the required L - tracking specifications 
required of the original system. 

Consequently, we start by considering a nominal system output specified 
by (w ao(t), z ao (t), y ao (t)) satisfying the following equations for the exponen- 
tially weighted system. 

z ao (t) = A n z a 0 (t) + B z w no (t) -B z <t> ao (y ao (t),t) (10.29) 
y ao(t) = C z Z ao (t) 

which is a completely known hypothetical plant corresponding to (10.26). 
<£ao(y ao(t),t) is a nonlinear design function, which is obtained by replacing 
the uncertainty in 4> a ( y 0 (*)> 7> *) with t* 16 nominal values. Equation (10.29) 
has the equivalent operator form 

y ao (<) = C z rB z w no {t) - C z rB z 4> ao (t)y ao (t ) + y ao (0) (10.30) 

where (f> ao {t)y ao (t ) = </> ao (y ao(<),^) and y ao (0) is the initial condition of the 
nominal system. It is easy to show that the system given by Eqn. (10.30) 
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is Lie stable (see Remark 2). 

Then, a comparison of the actual system and the nominal system (10.30) 
leads to the error operator equation 

y«(0 -y«(t) = c z rB z {w n (t)-w no {t))-c z rB z (<p a (t)y a (t)-<frao(t)y*o(t)) 

(10.31) 

on assuming y o (0) = y O o(0), from which we obtain 
y «(*) = C z rB z (w n (t)-w no (<)) 

-c z rB,(0 a (t)y o (t) - <t> wo (t)y a o(t)) + Yao(t) 

= ' l J r y a (t). (10.32) 

Thus, from (10.30) and remark 1, it immediately follows that W is an 
operator mapping Lie into itself. If w n (t) belongs to Li , it is easy to 
see that the satisfaction of the circle criterion of lemma 2 guarantees Li 
stability, because A n is Hurwitz and belongs to the sector [a,(3\. 

Hence, the operator V also maps L 2 back into itself. If w n (t) belongs to 
Lie, instead of Li, it is worth noting that the truncated function [w„(t)]r 
is still in L 2 , where T is the truncation time. Thus, up to the truncation 
time T, the input w n (<) is in L 2 , and so is the output y Q (t). 

Note that Eqn. (10.32) is a fixed point equation with the mapping T r : 
jr 2 — + £ 2 . Hence, we can utilize fixed point theory to obtain a criterion for 
the tracking performance. To this end, we summarize the local contraction 
mapping theorem in the following lemma: 

Lemma 4: Let y ao € (Y,||.||), Y is a Banach Space, and let B yw be 

defined as a closed sphere in Y with center at y ao given by 

B yw = {y a |||ya - y«o|| < <*n}, where a„ > 0 


(10.33) 
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Let ^ : Y — ► Y be a linear map. If there exists a constant k, 0 < k < 1, 
such that 

(0 ||¥y«-*ylll < *l|y«-ylll> v y „, yl g B yw 

( a ) \\*yao-y a o\\<an(i-k). 

Then, it follows that 

(i) ^ maps By w into itself 

(ii) ^ has a unique fixed point y* G B yw such that ’U'y* = y* a ■ 

Proof: See Martin (1976). 

Lemma 4, together with the fixed point equation (10.32), and the sphere 
specification on input-output tracking, lead to the following main result for 
guaranteeing Loo ■ tracking in the sense of spheres for the closed-loop system 
( 10 . 8 ): 

Theorem 2 : Let w(t) G R T and y (<) G R p denote the exogenous input and 
the plant output vector of the closed-loop system and let the design matrices 
Ki,K c ,H and the nonlinear design function fi 0 (t) be such that the closed- 
loop system satisfies the conditions stipulated in theorem 1. Let W 0 (t) 
(= [r^(<) d T 0 {t)\ T , d 0 (t) = 0), z 0 (t) and y 0 (t) be the known nominal input- 
output combination satisfying the nominal tracking model of Eqn. (10.30), 
and w(t)(= [r T (t) d r (t)] T ), z(<),y(<) be any combination satisfying the 
actual closed-loop equations (10.28). Then, for any reference input r(t) in 
the prespecified input sphere 

Q(r ,r 0 ,p Ti ) = {r (t) G I&| ||r(t) - r 0 (t)\\ Loo < M (10.34) 

and a bounded disturbance with norm ||d(t)||oo, there exists a unique plant 


CHAPTER 10. PROBLEM FORMULATION 


10-19 


output y(t) in the prespecified output sphere 

y,y.,&) = {yM € IU ||y(<) - y.Wlli™ < A) 

provided that 

a < @1 

Pl - &.+||d(0||oo + &ft + ft 

where 


(10.35) 

(10.36) 


A = ||C 2 (.sI — A„) -1 B^||// 00 , A n = al + A z , a >0 

/? 3 = sup ||[V y ^(y(t),7,t)]r||2 
7€9 

= sup ||(</>(y 0 (t),7,t) - </> 0 (yo(<),<))x||2 
769 

with 4> 0 {t) serving as the nominal value of 4>(t) at y 0 (t). 

Proof: First we consider the fixed point form (10.32) of the exponentially 

weighted system and obtain conditions for assuring that the mapping 'F is 
a local contraction over the exponentially weighted output sphere \\y a (t) - 
yao(<)l |2 < oin- 
Consider 

||(¥y a (0 - ^yl(0)r||2 = || - C z rB z (<t> a (t)y a (t) - <f> a (t)y' a (t )) T \\2 

< ||C 2 rB 2 ||, 2 ||(< 3 !> a (<)y 0 (t) - ^ a (<)yl(<))r ||2 

< fc||(ya(<)-yl(0)T||2 (10.37) 

where 


ll(6,(<)y.(*) -<k(*)yl(*))rlh 

k - ||C 2 rB 2 || i2 sup ||/ v (t) - v'«))tII 2 

y 0 #ya. ya,y a eB vw ,~fesscR r IfiyoUl y a \ l ))T\\2 
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= C,rB,|| a sup 


- 0.(t)yl(*))r||a 


ll(y.(0-y:(*))r||2 

From lemma 4, for $ to be a local contraction, we require 


(10.38) 


IK^yooCO - ¥y ao (<))r||2 

= ||C,rB,(w„(i) - w n0 (*))r - C z rB z (cf> a (t)y ao (t) - <f> ao (t)y ao (t)) T h 

< ||C,rB,||a{||(w n (<) - W no (t)) r || 2 + sup ||(</> a (t)ya O (0 - <£ao(<)yao(0)r||2} 

7G9 

< a„(l - lb). (10.39) 

Now, let ||(w n (t) — w no (f))r||2 < 0ni be a prespecified input sphere radius 
in L 2 and define 

01 = ||C z PB z ||,-2 = ||C z (sI — A„) _1 B z ||// 0O 

02 = SUp ||(^ o (0y«®(0 - <t>ao(t)yao{t))T\\2 

769 


03 = SUp 


ll(<^q(0ya(<) ~ ^a(<)yl(*))r||2 


yielding 

Expression (10.39) leads to 


Il(y«(0-ya(0)rll2 

k = 0i03 < 1- 


0\{0ni + 0?) < a n(l — k) 

which on combining with (10.41) yields 

0i < 


a T 


(10.40) 


(10.41) 


(10.42) 


(10.43) 


0ni + <*n03 + 02 

Hence from lemma 4, it now follows that inequality (10.43) is a suffi- 
cient condition for guaranteeing that ||(y Q (t) — y a o(f)):r||2 < <*„ whenever 
||(w n (f) - W no (<))r||2 < 0ni • 
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Note that inequality (10.43) assures L 2 tracking in the sense of spheres 
with respect to the input sphere radius P n i and the output sphere radius 
<* n . Since the design parameters are chosen to satisfy conditions of theorem 
1, we know that the original system is £«> - stable. Hence, we now use the 
L 2 - stability condition (10.43) to bound the Loo tracking error. 

Consider 

||(w n (<) - w«o(t))r||2 = (J Q |w„(0- w„ 0 (i)| 2 <ft) 1/2 

= (jf e 2 al | w (0 - w o (<)| 2 *) 1/2 

< ||w(t) - w 0 (t)|| 00 [^(e 2ar - 1)] 1/2 
= fl||w(<) - W 0 (t)||oo < Pni (10.44) 

where E = [i(e 2oT - 1)] 1/2 . Similarly, the relation between the variation 
of the output y a (t ) € L 2e in (10.10) and y 6 Loo in (10.8) is 

||(y.(t) - yao(<))r|| 2 < £||y(0 - y.(OII~ < • (10.45) 

Since we want the output of system (10.8) (measured by the Loo norm) to 
lie in the output sphere ||y(<) - y 0 (<)||oo < Po whenever the input w(<) is in 
the input sphere ||w(<) — w 0 (t)||oo < /?,, from (10.44) and (10.45), we obtain 

fini = Eft and (10.46) 

<*„ = Ep 0 . (10.47) 
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which is a sufficient condition for guaranteeing that 


|w (t) - w 0 (<)||oo < fa => \\y(t) - y 0 (f)||oo < A> • (10.49) 


Because w(<) in (10.49) consists of both the reference r(t) and the distur- 
bance d (f), by letting r 0 (t) and the null vector denote the nominal reference 
input and the nominal disturbance respectively, we compute 


|w(t) - W 0 (<)||oo = 


r (t) ~ r 0 (t) 
d(t) — 0 

< ||r(<) - r 0 (t)||oo + ||d(t)||oo < A (10-50) 


since the reference input r(t) is prespecified to be in the sphere ||r(f) — 
r 0 (t)||oo < pri, we let A = Pri + ||d(t)||<». Inequality (10.48) can then be 


rewritten as 

E(3o 

(h - E(Pr, + \\d(t)\\ 00 ) + Ep 0 p 3 + p2 
which is a sufficient condition for guaranteeing that 


(10.51) 


||r(f) - r 0 (t)||oo < Pri =>• ||y(<) - yo(0l|oo ^ Po (10.52) 


Remark 3: 

(i) In computing the upper bound on P\ specified by (10.51), a value 
for E is needed. E depends on the truncation time T. The system 
settling time, for instance, can be used for T. 

(ii) Pi is the norm of the linear operator C Z TB Z . Pi depends on 
the gains Ki, K c , H and the chosen positive number a. 

(iii) A, Pi and A can be computed as follows. 
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02 = SUp \\(<f> a (t)yao(t) ~ <t> ^o){t)T || 2 

7e9? 

< ||(e 0< )r|| 2 sup \\(<t>(y 0 (t),T,t) ~ &(y°(<),<))T h 

7G3 

= E sup ||(<^(y 0 (<), 7, <) - (f> 0 (y 0 (t ), 0)r||2 

7G3 

= (10.53) 


03 


ll(0a(Oya(*) ~ 

s - p ll(y.(<)-y:(*))r||2 

ll(e ot [<Ky(*)>7,*) -<?Hy'(0>~M)])rH2 

S “ P l|[e"(y(i)-y'(f)))T||2 

SUp||[Vy0] T ||2- 

* 


(10.54) 


Since </> € [a,/3], V 7 € $ and the bounds a,/ 3 are known apriori, we can 
simplify the computation of the inequality (10.51) by letting (f> 0 (y 0 (t),t) — 
||| a + (3\\y 0 yielding 

A = (H^J 1 )ll(y-(*»rlh - 0°' 88 ) 

Thus, the sufficient condition guaranteeing 


||r(<) - r 0 (<)||oo < 0rx =*• ||y(0 - yo(*)||oo < 0o (10.56) 


can be expressed as 

A < Si . (10.57) 

Pl ~ 0 rt + \\d(t)\\ oo + 0 o 0 3 + 0 4 K ' 

Remark 4: Although the parameters a and T do not explicitly appear in 
the inequality (10.57), the selection of these two parameters does affect it 
implicitly, since 0\ depends on a > 0 and 03,0s depends on T. To obtain 
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less conservative results, a and T can be chosen as small as possible. The 
smaller the values of a and T are, the smaller are the value of B\ and that 
of fh,04 in the inequality (10.57). 

10.5 Computation of the Gains K and H 

The sufficiency criterion (60) of theorem 2, essentially separates the uncer- 
tain effects from the nominal linear portion L = C Z TH Z ( = ||L||Woo)- 
Thus, the design problem now reduces to obtaining the state feedback gain 
matrices K = [Ki,K c ] and the observer gain matrix H such that /9j is 
close to its threshold value. To select K and H, we consider the follow- 
ing hypothetical linear system representing the transfer function matrix 
C z (sl - A n ) _1 B*. That is: 

iiii) = A„z ; (t) + B z w(t) (10.58) 

y(0 = c z z,(t) 


where 



al + A 

BK C 

BKj 


0 

G 

A n = 

— B c C 

al ■+• Ac 

0 

, B z = 

B c 0 


HC 

BK C 

al + A + BKi — HC 


L 0 

0 J 


C z = [C 0 0], 

By letting 

x,(f) 
x c ;(<) , 

x,(f) 


z i(t) = 
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Eqn. (10.58) can be expressed as the combination of state dynamics : 


X/(t) 

Xci(t) 


ol -(■ A 0 
— B c C al + A 

0 G 
B r 0 


X/(t) 

Xc/(<) 


+ 


B 

0 


u i(t) 


+ 


w(t) 


y(0 = [c o] 


x,(i) 

x c /(<) 


(10.59) 

(10.60) 


and of the full order observer for estimating the state x/(t): 


X/(0 

m 


(al + A)xi(<) + Bui(<) + H(y(<) - y(<)) 


Cx/(t) , u i(t) = [Ki K c ] 


x,(t) 

x c /(t) 


(10.61) 


In order to adjust the -norm of the hypothetical, linear system 
C z (sl — A n ) -1 B Z by proper choice of Ki, K c , and H we can use any 
^/oo -minimization algorithm available in the literature [Doyle et al (1989), 
Grimble (1988), Postlethwaite and Young (1988)] . The specific algorithm 
used in the present work was developed by the authors of the present paper, 
details of which can be found in Hwang et al (1990). 


Chapter 11 

An Algorithm for Jfoo 
Optimality 


The objective of the algorithm given below is the determination of the con- 
troller and the observer gain matrices so that the closed loop system remains 
stable and ||T W y Wh^ < v, where u is the threshold value of the if oo -norm 
of the transfer function matrix T W y from the hypothetical exogenous input 
w (t) to the hypothetical output y(<) of Eqn. (61). For a specific weighting 
p on the hypothetical control inputs u /(<), the if^-optimal control gains 
are determined by invoking theorem 3 and the existence conditions given 
below. 

Theorem 3: The ifoo-optimal control law for the weighted hypothetical 

y (*) 


controlled output vector 


pu,(t) 

[K, K c ] = 


of system (61) is given by 


ki 


(l\ 

B 


0 
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and H = — (I — hooki) -1 hoo[C 0] 


where ki is the positive definite solution of the algebraic Riccati equation 

iT 

ki + ki 


al -(- A 0 
— Bed al *f- Ac 


al "f A 0 
— BjC al -f Ac 


+ ki 




1 


1 

T 

1 



r\ 

0 

G 


0 G 

-(-) 

B 


B 



Be 

0 


1 

o 

y 

PQ 

i 

Uv 

0 


0 

) 


ki 


+ 


i T r 


C 0 


C 0 


= 0 ( 11 . 1 ) 
and hoo is the positive definite solution of the following Lyapunov equation 
[al + A]hoo + hoo[aI + A] T + GG T = 0 . (11.2) 


Proof: The proof is omitted since it appears in Hwang et al (1990). 

The necessary and sufficient conditions for the existence of the H 0 c - 
optimal control law given in theorem 3, such that ||T wy < u, are: 


(i) 


al 4- A 0 



0 

G 


0 

G 

T 

B 


B 

r\ 

— BcC al T Ac 



Be 

0 


Be 

0 


0 

rj 


0 

) 


i r 


■ 

T 

■ 

C 0 


C 0 


al + A 0 
— BcC al -(■ Ac 


(A + aI) T 0 
— GG t — (A -f- al) 


G dom(Ric) 


and 
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(ii) 


(Hi) 


ki 0 and hoo ^ ^ 


^mai (kxh 

oo ) ^ ^ 


where A mai (-) is the maximum eigenvalue of the matrix (•). 


A complete treatment of the above existence conditions can be found in 


Doyle et al (1988). 

The algorithm based on the existence conditions and theorem 3 can now 
be stated as follows for determining the relevant gain matrices. Henceforth, 
it is assumed that the design objective is H^TwyH//^ < 1 or equivalently 
IIT^Ik. < 1, for insuring the existence of the Jfoo-optimal control law. 
This transformation can be accomplished by scaling w(t) and/or y(t) such 
that the iZoo-norm of T wy is always less than or equal to 1. This scaling of 
the disturbance and controlled output terms requires the following trans- 
formation of the matrices involved in the hypothetical system dynamics: 


1 0 G 

-1/2 
— > V ' 

0 G 


B 

— ► i/ 1 ^ 2 

B 

o 

CQ 


B c 0 


0 


0 


,C -> n- 1/2 C. 

(11.3) 


We also define: 

IITwylk. £ 6 , (11.4) 

to facilitate the description of the algorithm below: 

Step 1: Select the (control) weight p and an initial value for v. The initial 
value of v is set at 1 1 T W y 1 1 , though alternate choices are acceptable 
[see Hwang et al (1990)]. 
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Step 2: Apply the aforementioned scaling and solve the Riccati equation 
(65) and the Lyapunov equation (66) for the scaled system. 

Step 3: Evaluate the gains ki and hoo corresponding to the unsealed sys- 
tem and calculate ||T W y||ffoo» i.e. the value of 6. 

Step 4: If ki > 0 and h<x, > 0, then go to Step 5, or else, 

if 6 > v, then decrease p and repeat Steps 2 through 4, 
or else, 

set i/ = <5 and repeat Steps 2 through 4. 

Step 5: If A mar (kihoo) < 1, then go to Step 6, or else, 

if 8 > v, then decrease p and repeat Steps 2 through 4, 
or else, 

set v = 6 and repeat Steps 2 through 4. 

Step 6: If 6 < v, then set v = 8 and proceed to Step 7, or else, 

if p < e , e el small number, then set v = 8 and go to Step 8, 
or else, 

decrease p and repeat Steps 2 through 6. 

Step 7: If u is less than the minimum value of v evaluated so far then 
repeat Steps 2 through 7 with a smaller value of v, or else proceed to 
Step 8. 

Step 8: Form the Roo -optimal controller gains Ki, Kcand the observer gain H. 
Moreover, ||T W y \\h^ ~ v. 
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Remark 5: The above algorithm is organized in such a way that it is 
versatile enough to accommodate other design specifications in addition to 
minimizing the i^oo-norm of T W y- 


Chapter 12 

Nonlinear Synthesis Procedure 


Based on the iT^-optimal control law given in the previous section, the 
controller gain matrices Ki,K c and the observer gain matrix H, needed 
for synthesis can be obtained for Lure type nonlinear systems through a 
step by step design procedure proposed below. This leads to a systematic 
placement of eigenvalues for the QPP problem for guaranteeing Too tracking 
in the sense of spheres for Lure type systems. 


12.1 Design Procedure 

First compute the threshold 

m = Wi + pwijT + Aft + fit (m) 

from the given specifications /3 0 , ||d(<)||oo? the nonlinear function 4>{t) 

and assumed initial guesses for a and T. m gives an initial upper bound 

°n I|L||h„ = ||C,(.I - A n )-'B,||„„. 
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Then use the Hoc - algorithm given previously in section IV to select 
a set of gain matrices K = [Ki K c ] and H so that ||IL» (1//^ < m. A new 
value of T can also be found by adopting the settling time of the closed- 
loop system. With these values of K = [Ki K c ] and H, we can determine 
a > 0 to satisfy the conditions of the circle criterion. With this newly 
found value of a and T, we can recompute the updated value mi of m. 
If the updated value mi « /?i, then the design is complete. If not the 
procedure is repeated with the newly found values of a and T. We have 
found the following step by step procedure to be effective for executing the 
design methodology developed in chapter 10. 

Step 1: Select values for a and T, and compute the threshold m. 

Step 2: Select p and use the H search algorithm 

until 1 = ||L|| Woo < m and determine K (= [Ki, K c ]) and H. 

Step 3: Form the closed-loop system. Adapt T to the system settling time 
and repeat step 1 to update m until T is close to the system settling 
time. 

If j3i > m, decrease p and repeat step 2. 

Step 4: If (jJ I * H.0' u>) satisfies the circle criterion, go to step 5. 

otherwise adjust a and go to step 1. 

Step 5: Repeat steps 2-4 

until (3\ is close to but lower than the threshold m. Then form the 
gain matrices K and H. 
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Two examples that utilize the above step by step procedure for executing 
the design method are given in the next chapters. 



Chapter 13 

A Simple Example 


Consider the uncertain, nonlinear system: 

x(i) = Ax(t) + Bu(t) + Gd(t) -f f(y(i), Am, t) (13.1) 

y(t) = C x(t) 

with 


A 




0 1 
0 -1 


B = 


0 

1 


, C = [l 0], 


G = 


0 

0.1 


0 

1 


{2 y(i) + (1.5 ± Am)y(f)smt} , 


where Am € [—0.05, 0.05] is the uncertainty, x(<) € R 2 are the states, y(t) 
is the output and d(t) is the external disturbance. 
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13.1 Design Objective and Performance Re- 
quirements 

The output of the above system is required to track a nominal sinusoidal 
input r 0 (t ) = sin 5 1 within a sphere of radius 0.2 for any reference input 
r(t) in a sphere of radius 0.1 centered at the nominal input, i.e., 

\\r(t) - r 0 (t)||oo < 0.1 =*► || y(t) - y o (0lloo < 0.2 (13.2) 

in the presence of the disturbance d(t) and the plant uncertainty f(y(t), Am, t). 


13.2 Design Execution 

To assure good steady state tracking, a servo compensator S c , denoted by 
(A c ,B c ,C c ), can be chosen as : 

x c (t) = A c x c (t) + B c e(t) 
y c (0 = C c x c (t) 


with the tracking error e(t) = (r(t) — y(t)), and |sl — A c | — s 2 + 25 — 0, 
which is the internal model of the signal sin 5t, yielding 


A c 


0 1 

-25 0 


B c = 


0 

1 


, Cc = [1 0] . 


(13.3) 
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With S c chosen as above, the following augmented nonlinear plant is ob- 
tained by combining S c and the plant model (13.1). 


X a (t) = 


0 

1 

0 

0 


0 

0 

-1 

0 

0 

X a (<) + 

1 

0 

0 

0 

1 

0 

-1 

0 

-25 

0 


0 


u(t) + 


0 0 

0 0.1 
0 0 

1 0 


w (<) + fi (y(t), Am,t) 


= A a x a (f ) + B a u(t) + G a w(i) + fi (y(<)» Xm,t) , 

y(t) = [1 0 0 0]xa(<), 


where 
X a (t) = 


x(t) 

Xc(<) 


, fi(y,Am,t) = 


f(y(t),Am,t) 

w (*) = 

r(t) 

0 

d(t) 



(13 


• 4 ) 


13.3 Computation of Gains K and H 


Invoking the ifoo-minimization algorithm described in Part I with p = 
0.0001, a = 10 — 8 , and 


z(t) = 


pui(t ) 
W ,e(t) 


, W s = S c , U,(<) = Kx a (t) 


(13.5) 


the controller gain matrices Ki, K c of the nonlinear controller, specified by 
Eqn. (9) of Part I and the observer gain matrix H of the nonlinear observer 
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specified by Eqn. (8) of Part I, are found to be 

Ki = [-299.9 -23.5] 

K c = [2506.2 1937.1] 

-0.0955 

H = 

-0.0046 

With the above observer gain matrix H and the controller gain matrix K, 
the closed-loop system becomes 

z (<) = A 2 z(t) + B z w(t) + B z </>(y(t), Am, t) (13.6) 

y(t) = C z z {t) , 


where 


A z 


B z 


C z 


0 . 

-23.5 

0 . 

0 . 

1. 

-24.5 
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0 . 

10 . 


(Am y(<) sint ) . 


It can be easily verified that the nonlinear element <j>(y(t), Am, t) belongs 
to the sector [-0.5 0.5]. 

By selecting T=0.5sec (the settling time), the Nyquist diagram u> 
H a (ju) satisfies the circle criterion, as shown in Figure 13.1, confirming 
that the closed-loop is L <*> stable. 

With the above choice of gains K = [K x , K c ] and H, the right hand side 
of the inequality (39) of Part I becomes 


m 


go 

Pri + ||d(0||oo + PoP3 + ^-WiVo^h 

0.2 

0.1 + 0.2 * 0.02727 + 0.05 * 0.363 “ 


1.61807 


and the computed value of f3 1 is 1.61454, which is clearly less than the 
threshold m. The resulting nonlinear control law is u(t) = K x x(t) + 
K c x c (f) - / 2 o(y(f),0, f2o(y(t),t) = 2 y(t) + 1.5 y(t) sin t. The above 

control law is optimal in the sense that the value of P\ is very close to 
its upper bound. To confirm the validity of the theorem, the closed-loop 
system of Eqn. (13.6) is simulated with the following results. 


13.4 Computer Simulation Results: 

The actual reference r{i ) is assumed to vary from the nominal reference 
input T 0 (t) with an amplitude of v 0 (t) — 0.1 for a period of 2 seconds (from 
0. seconds to 2.0 seconds). The histories of r 0 (t ) and r(i) are shown in 
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IM 



Figure 13.1: The Nyquist Diagram of H a (jio) and the Critical Disc D(a, (3). 
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Figure 13.2. Thus, || r(t) - r 0 (<)||cc = 0.1. Let Am = 0.05, and d(t) = 0. 
The nominal reference r 0 (t), and the actual output y(t) are shown in Figure 
13.3. As expected with the proposed nonlinear control law, as expected, 
the closed-loop system has low sensitivity to both the external disturbance 
and uncertainties. This is confirmed by Figure 13.4, showing that the 
actual output y(t) is close to the nominal one y 0 (t). Figure 13.5 shows that 
|| y - yolloo is less than 0.094. Clearly, good tracking can be achieved by 
using the proposed control law, u(t), whose control histories are shown in 
Figure 13.6. 

To check the validity of the design criterion another reference, r(t) = 
r 0 (<) -f 0.1 sinlOt, was simulated. The histories of r 0 (t) and r(t) are shown 
in Figure 13.7 with Am = 0.05, and d(t) = 0. The nominal reference r 0 (t), 
and the actual output y(t) are shown in Figure 13.8. Figure 13.9 shows that 
the tracking error ||y(t) — yo(0ll°° ^ ess 0-16, which is in the desired 

output sphere of ||y(t) - y 0 (*)||oo < 0.2. The resulting control history is 
shown in Figure 13.10. 
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Figure 13.5: Tracking Errors vs Time t. 
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igure 13.6: History of Control Input u(t). 
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Figure 13.7: Reference r(t) and Nominal Reference r 0 (t). 


Output 
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Figure 13.8: Output y(t) and References r 0 {t ) vs Time t. 




tracking error 
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V 


Figure 13.10: History of Control Input u(t). 



Chapter 14 

Two-Degrees-of- Freedom 
Robotic Manipulator 


Now, we consider the realistic problem of a two-degrees-of-freedom manip- 
ulator studied in Corless et al (1984). This arm can extend or retract and 
rotate about an axis as shown in Figure 14.1, where s is the distance be- 
tween the point load M and the pivot point 0. (xi(t), x 2 (t)) denotes the 
position of the mass center C of the arm in polar coordinates. l / i(t),l / 2 (t) 
are the force and the torque required to ensure that the point load M tracks 
a prescribed desired motion. (Vxi(0> ^ 42 ^)) denotes the actual position of 
the load M in cartesian coordinates, n and M are the arm and load mass 
respectively. J\ is the mass moment of inertia of the mechanism about the 
axis through 0; J 2 is the mass moment of inertia of the arm about the axis 
through C. a is the distance from C to the load M. The dynamic equations 
governing the motion of the manipulator can be written as 

(/i + M)x i(i) - R(xi(<),M)x^(<) = W) (14.1) 
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I(xi(t), M)x 2 (t) + 2R(x 1 (t),M)x 1 (t)x 2 (t) = V 2 (t) 

where 

R(xi(t),M) = pxi(t) + M( Xl (t) + a) 

I(xi(t),M) = J 1 + J 2 + px 2 1 (t) + M(x 1 {t) + a) 2 . 

Assume that all the parameters are known exactly except the payload 
mass M, which is subjected to known bounds, 0 < M < M < M, denoted 
as M G [M . M]. Let M° be the nominal payload mass (= + M)). 

By defining 

til (0 = (/i + 

ua(0 = (/(*i(t),M))- 1 V 2 (t) 

as the new input variables, and 

/.I (0 = + 

U(t) = -(,I{x I (t),M))-'2R(x,(t).M)x 3 (t)x,(t) 

as the modified nonlinear functions with 




Xi(t ) 


0 

Ul(t) 

U2 (t) _ 

, x(t) = 

X 2 (t) 

Xz(t) 

, and f(x(t),M,t) = 

0 



X 4 (t) _ 


_ fn2(t) _ 


the system (14.1) can be rewritten in the following state equation form: 

x(t) = Ax(t) + Bu(<) + Gd(t) + f(x(<), M, t) (14.2) 

y (0 = Cx(t) 


CHAPTER 14. 2 DOF Robotic Manipulator 


14-3 


Va 2 



Figure 14.1: Robot Manipulator of Two-degree-of Freedom. 
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with 


A 


C 


0 0 10 
0 0 0 1 
0 0 0 0 
0 0 0 0 

10 0 0 
0 10 0 


0 0 
0 0 
1 0 
0 1 


G = 


0 0 

0 0 

0.1 0 

0 0.1 


where u(t), x(t), y(t), d(t) and f(.) are respectively the control, the state, 
the output, the disturbance and the nonlinear uncertain element of the 


system of Eqn. (14.2). 

Parameter values used for computations are M = 45 kg, M = 55 kg, 
Jj = J 2 = 100 kg.m 2 and a = 1 m, with uncertain mass M £ [45 kg, 55 kg]. 


14.1 Design Objective and Performance Re- 
quirements 

Suppose the payload M is required to track the semi-circle described by 
[Y A1 (t),Y A2 (t)) = [0 m, 1.5 m] , <<0 

[Y Al {t),Y A2 (t)] = [(*!(<) + a)cos x 2 (t) m,(x!(t) + a)sin x 2 (t) m] 

0 < t < 3.1416 sec. 

[Y AX {t),Y A2 {t)\ = [0m.-l.5m], t > 3.1416 sec. 



CHAPTER 14. 2 DOF Robotic Manipulator 


14-5 


In terms of polar coordinates the above desired motion, becomes 

MO, yr 2 ( 0 ] = [°- 5 m ’ 0 ra< ^ » * - 0 sec - 
[y rl (<),y r2 (0] = [0.5 m, t rad] ,0 <<< 3.1416 sec. 

MO.MOl = [0.5 m, 3.1416 rad\ ,t > 3.1416 sec. 
where y rl (t), y r 2 (0 represent the desired motion of the system output y(t)(= 

[yi(t),y 2 ( 01 r ) in E( i n - ( 14 - 2 )- 

The performance requirement is that 

||r(0 - r 0 (<)||oo = 0 =» ||y,(0 - y o ,(0ll°° < 0.14, * = 1,2 , (14.3) 

where the nominal reference input r 0 (<) = [y r i(t),yr2(t)} T - 

14.2 Design Execution 

For steady-state tracking of the desired motion, a servo-compensator S c 
denoted by (A c , B c ,C c ) is chosen as : 

x c (t) = A c x c (i) + B c e(f) 
y c (t) = C c x c (t) , 


where the tracking error e(f) = (r(<) - y(<)), 
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The above servo- compensator S c is chosen to realize the internal model, 
which is specified by the desired motion, so that the desired trajectories 
can be followed. With S c chosen as above, the augmented nonlinear system 
becomes 



Now we invoke the design procedure outlined previously in Part I. 


14.3 Computation of the gains I< and H 

Choosing 

x(t)= PU ' W W s = S c , e(t) — r(t) — y(t) (14.6) 
_ W ,e(t) _ 

where u i(t) = Kx a (<), with p = 0.001 and a = 10 -9 , the state feedback 
gain matrix K = [Ki K c ] and the observer gain matrix H are computed, 
yielding : 

-202.9 0. - 20.2 0. 1020.4 0. 0. 

0. -108.4 0. -14.7 0. 1005.9 466.9 
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-0.4571 

0.0000 

0.0000 

-0.4501 

-0.1022 

0.0000 

0.0000 

-0.1007 


14.4 Formulation of the Closed-loop System 


The closed-loop system can now be expressed as 

i(t ) = A z z(t) + B z w(t) + B z <f>(x(t),M,t) (14.7) 

y (t) = C z z(t) , 


where A 2 ,B 2 , and C z are shown in Table 14.1, and 

<f>(x(t),M,t) = 

0 0 
0 0 

'♦O 

£ 

0 

1 

-to 

10 0 
0 10 



fi(x, M,t) = 

_ _ 


-(I{x 1 (t),M))- 1 2R(x 1 (t),M)x 3 (t)x 4 (t) _ 

t— * 

o 

III 

fiol (^) 


(p + M^-'Rix^M^xKt) 

_ f2o2{t) _ 


-(I(x i (t),M°))- 1 2R(x 1 (t),M°)x 3 (t)x 4 (t) 


with the nominal mass M° = 50 kg. 

Remark: Note that f lo (x,<) is the well known computed torque. 


With a = 10 -9 and T = 3.1416 sec, Figures 14.2 and 14.3 show that 
the Nyquist locus u> t-+ H a (ju>) satisfies the MIMO circle criterion, -proving 
that the closed-loop system is Loo stable. 
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erion (for Sector = [«i, /?i])- 
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IM 



Figure 14.3: Satisfaction of the Circle Criterion (for Sector = [a 2 , fa])- 
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The numerical values of the parameters in inequality (39) of Part I are 
computed as follows: 


/?3 


A 


sup ||[Vy<£(y(f),M,t)]r||2 
A/€$ 

r d<£n dcf>2\ i || 
rs — “5 — It 2 

oyi oyi 

0.018 


sup \\(4>(y 0 (t),M,t) - {<t> 0 (y 0 (t),i))T \\2 - 0.056 . 

Mea 

(14.8) 
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TABLE 14.1 The Closed-Loop System Matrices of the Robot 

Manipulator 


A z (11 x 11 matrix) = 


10 3 * 


Columns 1 thru 8 


0.0000 

0.0000 

0.0010 

0.0000 

0.0000 

0.0000 

0.0000 

0.0010 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0010 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0010 

0.0000 

0.0000 

0.0005 

0.0000 

0.0000 

0.0000 

0.0000 

0.0005 

0.0000 

0.0000 

0.0001 

0.0000 

0.0000 

0.0000 

0.0000 

0.0001 

0.0000 

0.0000 

Columns 9 thru 11 


0.0000 

0.0000 

0.0000 


0.0000 

0.0000 

0.0000 


0.0000 

-0.0202 

0.0000 


-0.1084 

0.0000 

-0.0147 


0.0000 

0.0000 

0.0000 



0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0204 

0.0000 

0.0000 

-0.2029 

0.0000 

1.0059 

0.4669 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0010 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0005 

0.0000 

0.0000 

0.0000 

0.0000 

1.0204 

0.0000 

0.0000 

-0.2030 

0.0000 

1.0059 

0.4669 

0.0000 
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TABLE 14.1 (cont.) 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0010 

0.0000 

-0.0005 

0.0000 

0.0010 

0.0000 

-0.0202 

0.0000 

-0.1085 

0.0000 

-0.0147 


B z (11 x 

4 matrix) = 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.1000 

0.0000 

0.0000 

0.0000 

0.0000 

0.1000 

1.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

-0.4571 

0.0000 

0.0000 

0.0000 

0.0000 

-0.4501 

0.0000 

0.0000 

-0.1022 

0.0000 

0.0000 

0.0000 

0.0000 

-0.1007 

0.0000 

0.0000 
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TABLE 14.1 (cont.) 


C z (2 x 11 matrix ) = 

Columns 1 thru 8 

1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

Columns 9 thru 11 

0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 


Thus, the threshold m becomes 

go 

m ~ Ai + ||d(0IU + A& + A 

9 3923 

_ 0. + 0.14 * 0.018 + 0.056 

and the computed value of ft is 2.34285, which clearly satisfies the in- 
equality (39) of Part I. The resulting nonlinear control law is u(t) = 
K 1 x(f) + K c x c (t) — fi 0 (t). The above control law is optimal in the sense that 
the value of ft is close to its upper bound. Prom theorem 2 of Part I, the 
closed-loop system is now guaranteed to have the required tracking perfor- 
mance. Effectiveness of the above controller is confirmed by the following 
transient response simulation results. 
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14.5 Computer Simulation Results: 

Let d(t) be an impulse disturbance. With the lower bound mass M = 45 kg, 
the plant outputs y (<) are shown in Figure 14.4, which shows that the plant 
outputs y\(t) and j/ 2 (t) closely track the reference inputs T\ (t) = 0.5 m and 
— t radians respectively. The maximum tracking errors between the 
plant outputs and the nominal outputs are 0.006 m and 0,0015 radians. 
The errors are relatively small in comparison to the imposed output sphere 
radius jd 0 = 0.14. The superiority of the proposed control law for distur- 
bance rejection and its low sensitivity to uncertainties can clearly be seen 
from the trajectories of the masses M = 45 kg , M = 50 kg, and 

M = 55 kg in the Cartesian coordinates, which are shown in Figure 14.5. 
The histories of the required control force Vi(<) and the torque V 2 (t) for 
M = 45 kg are shown in Figure 14.6 along with those for the nomi- 

nal mass M = 50 kg. Due to the desired high-speed tracking, it is not 
surprising that the required control inputs are relatively high. 



y^t) (m) Mid y a (£) ( radians ) 
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Figure 14.4: The Plant Output y(t) vs Time t. 
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)• 




V t (e) (N) and V,(t) (Nm) 
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Time (sec.) 


Figure 14.6: Histories of Control Force Vi(t) and Torque V 2 (t). 







Chapter 15 

Summary and Conclusions 


In this paper the problem of tracking in the sense of spheres, for a class 
of multivariable nonlinear uncertain systems of the Lure type, is consid- 
ered. Design criteria developed, extend the circle criterion developed in 
Jayasuriya & Kee (1988) for the L 2 - problem to L 0 0 - tracking. Controller 
synthesis is reduced to a quantitative pole placement scheme which is sub- 
sequently solved using the H <*>■ -formalism. The Hoo -formalism provides a 
framework for incorporating bandwidth issues, so crucial in any physical 
implementation of any resulting control law. However, in light of the op- 
erator methods used in the problem formulation, resulting designs can be 
somewhat conservative. This, is expected since any Loo -function within the 

S 

given input sphere is admissible. 

Future work will consider alternate approaches for reducing the conser- 
vativeness of the proposed design criterion. Also, it would be very useful to 
have a complete characterization of the worst input lying in a given input 
sphere for the synthesized closed-loop system. The worst input here is the 


15-1 



CHAPTER 15. SUMMARY AND CONCLUSIONS 


15-2 


input which leads to the maximum error in the actual output as measured 
with respect to the nominal output. Criteria for synthesizing controllers 
for more general nonlinear system are also needed. For example, it would 
be appropriate to study polynomial systems, discontinuous systems and 
systems in which the control does not enter the state equation in a linear 
way. 

Two examples for the synthesis methodology presented in Part I of this 
paper are given. The simulation results confirm the effectiveness of the 
proposed design methodology. The -framework used in the synthesis 
procedure can be adapted to incorporate bandwidth limitations. The lat- 
ter can be accomplished by adjusting the weighting p and the filter matrix 
W a needed in the Hoc- algorithm. It should be emphasized, however, that 
incorporating such bandwidth limits may result in having to relax the track- 
ing error specifications. The relationship between the tracking error and 
the bandwidth is implicitly given in the design criterion (60) developed in 
Part I of the paper. Qualitatively, a large bandwidth would yield a small 
.Hoc -norm and vice-versa. Moreover, for a given bandwidth limitation, a 
controller can be found, whenever the minimum ffoo-norm with an appro- 
priate set of p and W a , is smaller than the threshold value m (see Eqn. 
(60) of Part I). If this is not the case the same inequality can be used to 
ascertain a new output sphere size. 
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